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Abstract 

The Bogoliubov-Dirac-Fock (BDF) model is the mean-field approxi- 
mation of no-photon Quantum Electrodynamics. The present paper is 
devoted to the study of the minimization of the BDF energy functional 
under a charge constraint. An associated minimizer, if it exists, will usu- 
ally represent the ground state of a system of A'^ electrons interacting with 
the Dirac sea, in an external electrostatic field generated by one or several 
fixed nuclei. We prove that such a minimizer exists when a binding (HVZ- 
type) condition holds. We also derive, study and interpret the equation 
satisfied by such a minimizer. 

Finally, we provide two regimes in which the binding condition is ful- 
filled, obtaining the existence of a minimizer in these cases. The first 
is the weak coupling regime for which the coupling constant a is small 
whereas aZ and the particle number A'^ are fixed. The second is the 
non-relativistic regime in which the speed of fight tends to infinity (or 
equivalently a tends to zero) and Z, N are fixed. We also prove that 
the electronic solution converges in the non-relativistic limit towards a 
Hartree-Fock ground state. 



1 Introduction 

The relativistic quantum theory of electrons and positrons is based on the free 
Dirac operator [14j . which is defined by 

3 

D° = -ihc ^akdk+ mc^j3 -.^ -ihca-W + mc^ (3 (1) 

k=l 

where a. = (ai,a2,a3) and (3 — ( ^ '^r ) ' '^^ I ^ '^k 
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"' = (^1 o)' o')' ^'^(o -i)- 

We follow here mainly the notation of Thaller's book :5Sj. In (IT|), ?i is Planck's 
constant, c is the speed of light and m is the mass of a free electron. For the sake 
of simplicity, we shall use in the following a system of units such that h = m = 1 . 
Unless otherwise specified, we shall also assume that c — 1, in which case an 
additional parameter will appear in front of the interaction potentials, a = e^, 
where — e is the bare charge of a free electron. 

The operator acts on 4-spinors, i.e. functions ip S L^(R^,C^). It is self- 
adjoint on L2(k3^ (^4)^ domain H^{m.^, <C^) and form-domain TJ^/^^j^s^ ^4-,^ 
Moreover, it is defined to ensure {D^Y = — A + 1. The spectrum of is 
cr(Z)") = (— oo, —1] U [l,oo). In what follows, the projector associated with the 
negative (resp. positive) part of the spectrum of will be denoted by P° 
(resp. Pi): 

P° X(-co,o)P°), Pi X(o,+oo)P°). 

We then have D^P° = P°D° = -\D°\ P^ and D°Pl = PlD° = \D°\ Pi. 

Compared with the non-relativistic (Schrodinger) models in which —A/2 
appears instead of D", the main unusual feature of the relativistic theories is 
that (j{D'^) is not bounded from below. Indeed the free Dirac operator ([T]) was 
proposed by Dirac in 1928 [14] to describe the energy of a free relativistic spin- 
1/2 particle like an electron. In order to explain why negative energy electrons 
are never observed, Dirac made the assumption [Kl UHl [ITj that the vacuum is 
filled with infinitely many virtual electrons occupying all the negative energy 
states so that, due to the Pauli principle, a physical free electron cannot have a 
negative energy. This model is commonly called the Dirac sea. Mathematically, 
the free vacuum is identified with the projector P". 

With this interpretation, Dirac was able to conjecture the existence of the 
positron (the anti-electron, which has a positive charge), which is seen as a hole 
in the vacuum and was discovered in 1932 by Anderson ^j. He also predicted 
interesting new physical features as a consequence of his theory [TH Uni fT7] . 
which were experimentally confirmed later. First, the virtual electrons of the 
Dirac sea can feel an external field and they will react to this field accordingly, 
i.e. the vacuum will become polarized. From the experimental viewpoint, vac- 
uum polarization plays a rather small role for the calculation of the Lamb shift 
of hydrogen but it is important for high-Z atoms [44 and it is even a crucial 
physical effect for muonic atoms IHl [H]. Second, in the presence of strong 
external fields, the vacuum can acquire a nonzero charge, a phenomenon which 
is related to the spontaneous creation of electron-positron pairs [13 H?! HH 112] • 

On the other hand, many models which are commonly used to describe rela- 
tivistic particles do not take the vacuum polarization effects into account. This 
is for instance the case of the (mean-field) Dirac-Fock theory which is the rela- 
tivistic counterpart of the well-known Hartree-Fock model and was proposed by 
Swirles [55^. The Dirac-Fock model suffers from an important defect: the corre- 
sponding energy is not bounded from below, contrary to the Hartree-Fock case, 
and this leads to important computational difficulties (see [S] for a discussion 
and detailed references). From the mathematical viewpoint, one can prove that 
the Dirac-Fock functional has critical points which are solutions of the Dirac- 
Fock equations [121 SS] , but these critical points have an infinite Morse index, 
and the rigorous definition of a ground state is delicate [501 HI] ■ 
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It was proposed by Chaix and Iracane [101 that these difficuhies could be 
overcome by incorporating the vacuum polarization effects in the theory, i.e. 
by considering the coupled system 'Dirac sea + real electrons' instead of the 
electrons alone. Starting from Quantum Electrodynamics (QED) and neglecting 
photons, they derived a model called Bogoliubov-Dirac-Fock (BDF), in which 
the real particles are coupled to the Dirac sea. The main advantage of this 
theory is that the energy of the model is now bounded below, leading to a clear 
definition of the ground state. 

The Chaix-Iracane model was first mathematically studied in the free case 
by Chaix, Iracane and Lions in [TT] and then by Bach, Barbaroux, Helffer and 
Siedentop in [I] . The external field case was rigorously defined and studied by 
the authors of the present paper in [271 HH] ■ Chaix and Iracane derived their 
functional under Dirac's assumption that without external field the vacuum is 
given by P^.. This choice is not physically correct: it corresponds to neglecting 
the interaction between the virtual electrons. This deficiency was recently over- 
come by Hainzl, Lewin and Solovej [23] who used a thermodynamic limit applied 
to the QED Hamiltonian restricted to Hartree-Fock states (mean-field approxi- 
mation), in order to define the free vacuum. Doing so, they obtained a slightly 
different translation-invariant projector T-"" , solution of a certain self-consistent 
equation. Then, they showed by the same thermodynamic limit procedure that 
in the external field case the BDF model should better rely on this new free 
vacuum instead of Dirac's choice P° used by Chaix and Iracane. Note that the 
projector had been constructed earlier by Lieb and Siedentop ^j, but the 
existence proof and the physical interpretation were different. 

The Bogoliubov-Dirac-Fock model is a very promising theory: it is well- 
justified physically, it is better behaved than the usual Dirac-Fock model and 
it leads to new mathematical problems which are interesting in themselves. In 
particular, a state of the system always contains infinitely many particles (the 
real and the virtual ones). This property which raises serious mathematical 
difficulties is shared with other quantum models, to which a similar study could 
be applied. 

The purpose of the present paper is to continue the study which was started 
in [271 [251 [25] . In the BDF model, the state of the system is represented by an 
orthogonal projector of infinite rank 

i>l 

where {ifi)i>i is an orthonormal basis of Ran(P). The projector P should 
be seen as the one-body density matrix of the following formal wavefunction 
depending on infinitely many variables 

^ = tpi Atp2 A--- Aif, A-- - , (2) 

which is a kind of infinite Hartree-Fock state. Here A denotes the usual wedge 
product of functions in L'^{M.^, C*). The projector P represents the whole system 
consisting of both the real and the virtual particles of the Dirac sea, but there is 
no distinction between them a priori. It is only for the solution of the problem 
that the real particles will be identified and separated from the virtual ones. The 
BDF energy f is a nonlinear functional of the variable P, which is bounded 
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below on an appropriate set which will be defined later. The expression of f 
depends on a fixed external charge density ly. 

In [27ll28j . we proved the existence of a global minimizer of the BDF energy, 
interpreted as the polarized vacuum in the electrostatic field V = —av * | • 
created by the external density v. Here a = is the bare coupling constant 
and — e is the bare charge of an electron. The density v has no sign a priori but 
when V > one may think of it as the density of a system of nuclei. In this 
paper we study the minimization of under the constraint that the charge 
of our state P is equal to —eN where A'^ is some integer. Of course the total 
charge of a state of the form ^ is formally infinite since VE* represents infinitely 
many negatively charged particles. But one can define the difference between 
the charge of P and the (infinite) charge of the free vacuum. It is this difference 
which is fixed to —eN. A precise definition will be given below. 

This charge constrained minimization problem is much more delicate than 
the global minimization of [571 [5Sj. When > and iV > 0, a minimizer in the 
—eN charge sector will usually represent the state of N electrons interacting 
with the polarized vacuum and the external field V = —aiy* | • |^^- As usual for 
Hartree-Fock type theories [HI [351 US] , one does not expect that this minimizer 
will always exist. Indeed, if v is not strong enough to bind the N electrons 
together with the polarized vacuum, there should be no minimizer. On the 
other hand, if v is too strong, some electron-positron pairs could be created. 

Let us denote by E"{N) the infimum of the BDF energy f in the charge 
sector —eN and in the presence of the external density v (a precise definition will 
be given in Section [2.3. ip . Our main result (Theorem [T|) will be the statement 
that all the minimizing sequences for E'^(N) are precompact if and only if a 
HVZ-type inequality holds: 

VifeZ\{0}, E^iN) < E^iN ~ K) + E'^{K). (3) 

Inequalities like ^ are very common in the study of linear [32 [Ml [S31 US] 
and nonlinear [121 [43] systems. Assume iV > for simplicity. Then E'^{N) 
is the infimum of the energy of a system of TV electrons coupled to the Dirac 
sea. When < K < N, ^ means that it is not favorable to let K electrons 
escape to infinity, while keeping N — K electrons near the nuclei. When K < 0, 
it means that it is not favorable to let \K\ positrons escape to infinity, while 
keeping + electrons near the nuclei. When K > N, it means that it is not 
favorable to let K electrons escape to infinity, while keeping K — N positrons 
near the nuclei. When a is small enough and > 0, it will be shown that the 
separation of electron-positron pairs is not energetically favorable, so that one 
just needs to check ^ ior K — 1, 2, N. 

From a mathematical point of view, proving the compactness of minimizing 
sequences assuming ([3]) is a subtle task. Indeed, the fact that our main variable 
P is a projector of infinite rank complicates a lot the study of minimizing se- 
quences, for instance compared to the Hartree-Fock case [HI [53^ in which only 
finitely many particles are described. In particular, it is not obvious at all to 
localize our state P in space in order to decouple the electrons staying close to 
the nuclei from those which escape to infinity. These complications are conse- 
quences of the vacuum effects. Similar issues are encountered in the study of 
other models describing systems of infinitely many particles. It is our hope that 
this work will provide a better understanding of these other models too. 
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When ([3]) holds, there exists a minimizcr in the charge sector —eN. This is 
an orthogonal projector P satisfying some nonlinear equation of the form 



P = X(-oo.,,] {D' 



\-\-' + Vp) 



(4) 



where Vp is an operator depending on P itself and /i is a Lagrange multiplier 
associated with the charge constraint, interpreted as a chemical potential. The 
same equation was obtained in [271 [2S HI! for the vacuum case (global mini- 
mization), but with fj, = 0. When the external density v is not too strong and 
> 0, then it will hold > and the operator 13° — aiy * \ ■ \ + Vp will have 
exactly N eigenvalues (counted with their multiplicity) in (0, /i]. In this case P 
can be written 



P = Xi-oofi]{D°-a'^*\-r' + Vp)+Xio,t.]{D°-a,^*\-r' + Vp) 



Formula (O allows to distinguish the "real" electrons (represented by the or- 
bitals {fn)n=i) from the self-consistent polarized vacuum Pvac- As explained 
in Section B.3.1[ the orbitals {fi, ...,ipn) are solutions of a Dirac-Fock type 
system of equations [T^, in which the mean- field operator is perturbed by the 
self-consistent vacuum polarization potentials. 

In the present work, we shall provide two regimes in which the condition 
^ holds, and therefore for which a BDF minimizer exists in the charge sector 
—eN. The first is the weak coupling regime in which the coupling constant 
a <C 1, but ai^ (hence aZ) and N are both fixed (Theorem [2]). The second is 
the nonrelativistic regime c S> 1 with v and N fixed (Theorem [3]) . In the latter 
case, we also prove that the N orbitals {(pi, (Pn) of ^ converge to a ground 
state of the nonrelativistic Hartree-Fock functional [HUl^ as c ^ oo. A similar 
result was already obtained by Esteban and Sere in the Dirac-Fock case |20| . 

The paper in organized as follows. In the first section, we define properly 
the model and state our main results. Section [3] is devoted to the proof of some 
preliminary results which will be needed throughout the paper. The last three 
sections are then devoted to the proofs of Theorems [l] [2] and [31 

Acknowledgement. The authors acknowledge support from the European 
Union's IHP network Analysis & Quantum HPRN-CT-2002-00277. M.L. and 
E.S. acknowledge support from the project "ACCQUAREL" NT05-4_44652 
funded by the French National Research Agency GIP-ANR. 

2 Model and main results 

2.1 The mean-field approximation in no-photon QED 

We start by recalling briefly the physical meaning of the model, as explained in 
[29j . As mentioned in the introduction, the state of our system is represented by 
an infinite-rank orthogonal projector P, which is seen as the density matrix of an 
infinite Hartree-Fock state ^ . We recall that although P should be interpreted 
as the state of the coupled system 'real particles -1- vacuum', there is no canonical 



N 




(5) 



n=l 
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distinction between the real and virtual particles a priori. For N electrons, this 
would correspond to a decomposition of the form P = Pvac + 7 where 7 is an 
orthogonal projector of rank N satisfying Pvac7 = 7-Pvac = (for N positrons, 
this becomes P — Pvac — 7)- There are infinitely many such decompositions for 
a given P and a given N. But for the solution of our equation, a particular 
decomposition may be chosen in a natural way. 

The energy of the system in the Hartree-Fock state P can be deduced from 
the QED Hamiltonian formalism in Coulomb gauge and neglecting photons, see 
[29j . The energy functional is formally 



V = —ai'* I ■ 1^^ being the external electrostatic potential created by the density 
of charge v. If > it can for instance represent a system of nuclei in a molecule, 
but in most of the results of this paper the sign of needs not be fixed. Also 
we shall not allow pointwise nuclei and v will essentially be an L^^^ function. 
It is a well known phenomenon in QED that pointwise nuclei create spurious 
divergences (see e.g., [35] ) . But the regularity assumption on v is not really a 
restriction from the point of view of physics: point-like nuclei do not exist in 
nature. In ([7]), pr is the density defined formally as prix) ~ tr^i (T (x , x)) , and 
a is the bare coupling constant, a = e^. 

In QED, a global minimizer of P 1-^ ^qed(^ ^ 1/^) represents the vacuum, 
whereas other types of states (for N electrons for instance) are obtained by 
minimizing this functional with a charge constraint. The charge is formally 
defined as 



The subtraction of half the identity in ^ and ([5]) is a kind of renormalization 
which was introduced by Heisenberg 31j and has been widely used by Schwinger 
(see [53, Eq. (1.14)], 54, Eq. (1.69)] and ^ Eq. (2.3)]) as a necessity for a 
covariant formulation of QED. 

Of course, the expression of the mean-field QED energy ^ is purely formal: 
if P is an orthogonal projector in infinite dimension, P — 1/2 is never compact 
and therefore £qej:,{P — 1/2) is not well defined. Even the density of charge 
P[p-i/2] is not a well defined object. For this reason, it was proposed in 23] to 
use a thermodynamic limit in order to give a rigorous meaning to the minimiza- 
tion of ([7]): the idea is to define the model in a bounded domain in space, and 
a cut-off A in Fourier space. This was done in in a box Cl = [—L/2; 
with periodic boundary conditions, and cutting the Fourier expansion outside a 
ball of radius A. Then, the minimization in Cl makes perfectly sense (i^(R'^, C^) 
has been replaced by a finite-dimensional space), and one can study the limit 
of the sequence of minimizers when L — > 00. 

In [29] . this technique was used to define properly the free vacuum and 
justify the validity of the BDF functional. Notice that the ultraviolet cut-off A 



^^fQED(^-l/2) 



(6) 



where 
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is fixed and will not be removed: it is well-known that QED contains problematic 
ultraviolet divergences which are difficult to deal with. We therefore introduce 
the following functional space 

i3A := {/ e L\R\,C^), supp(7) C 5(0, A)} . 

Notice that Sja is contained in the domain H^{M.^, C*) of and that D'^Sja ~ 
S^A- In the following, we still denote by D° its restriction to i^A. Taking v = in 
([7]) (free case) and studying the thermodynamic limit L oo, the free vacuum 
was obtained in 29J. It is a translation-invariant projector satisfying the 
Euler-Lagrange equation 

V" = X(_oo,o) {-D") , 
.o_^o iVl~l/2)ix,y) (9) 



\x - y\ 



The operator which appears in ([5]) is a translation-invariant operator taking 
the following special form [iQlESj, in the Fourier space, 

V'^{p) = a-Upg,i\p\) + goi\p\)P, ^ p/\p\. (10) 

Here gi and go are real and smooth functions satisfying 

X < gi{x) < xgo{x). (11) 

Note the self-consistent equation © was already solved by Lieb and Siedentop 
[40j . but their interpretation was not variational. They used a fixed point ap- 
proach valid when a log A < C. In [29j, the free vacuum solution of ([9]) is 
constructed as a minimizer of the energy per unit volume for any value of the 
ultraviolet cut-off A, and under the condition < a < A/ir. This last inequality 
is related to Kato's inequality \x\~^ < 7r/2|V|. Hence, in the whole paper we 
shall assume that < a < A/tt. We use the following notation = 1 — V'^. 

The mean-field approximation in no-photon QED is therefore very close 
to the original Dirac's picture of the free vacuum, the latter being described 
as an infinite rank spectral projector associated with the negative spectrum 
of a translation-invariant Dirac-type operator. However, it does not correspond 
exactly to the original ideas of Dirac: when a 7^ 0, 7^° is different from Even 
the free vacuum is solution of a complicated nonlinear equation This 
is because the interaction between the virtual particles is taken into account, 
similarly to the real ones. The Dirac picture is only recovered in the non- 
interacting case a = 0. 

It is important physically that the so-obtained free vacuum is invariant by 
translations. This means that the density of charge P['po_i/2] is (formally) 
constant. More precisely, the subtraction of half the identity allows to obtain a 
vanishing density, p^-po -1/2] = 0. By (jlOp . we have 

r° (v) - 1/2 = gi(bl)Q-^p + go(|p|)/3 
2x/5i(H)2+5o(H)2 ' 

from which we infer that trc4[(7'° - l/2)(p)] = for any p e -0(0, A), the 
Pauli matrices being trace-less. Thus the (constant) density of charge of the 
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free vacuum vanishes: 

P[po_i/2] = (2^)"' / trc4(7'"(p) - 1/2) dp = 0. 

JB{Q,A) 

This formally means that 

"tr{Vl-l/2)= f p[po_i/2](a:)dx = 0" (12) 

and therefore that the free vacuum is not charged. 
As a consequence of (ITUl) . the spectrum of I?" is 

a{V") = {±./^^iM)^TM^, P e B{0,A)} . 

It has a gap which is greater than the one of Z?", by (fTTj) : 

1 < to(q;) := mincrdP^I). (13) 

In Lemma Il9l below, we shall prove that when a <ti 1, then m{a) = go(0). We 
conjecture that this is true for any < a < 4/tt. Notice that the following 
expansion is known [101 US]: 5o(0) = 1 + f arcsinh(A) + O(a^). 

Once the free vacuum is defined, in the external field case v one can 
measure the energy of any state P with respect to the (infinite) energy of the 
free vacuum Pi . The Bogoliubov-Dirac-Fock energy is formally defined as 



■£:-(P_7,0) ^ £-^j,j,(P-1/2)-£^ed(^--1/2) (14) 

dx 

\Q{^,y)\' 



tr (P«g) + / F(.),g(.)dx + Uf '-^^^f^dxdy 
Jw 2 J Jk^xv? f - y\ 



a 



2 yyR3xR3 \x-y\ 



-dx dy " , 



where we have introduced Q = P — = [P — 1/2) — [V^ — 1/2) and used the 
definition of the self-consistent free Dirac operator , see (ITU)) . This functional 
has been mathematically defined and studied in [371 HH] j but with P° as refer- 
ence. In [21], it was proved that any sequence of global minimizers of the full 
QED energy in a finite box of size L converges (up to a subsequence) to a global 
minimizer of the BDF functional as L — > oo, justifying the formal derivation 
(fH|) . Notice that by , the total charge of our state is now formally given by 

"- etr(P- 1/2) = -etr(Q) = -etr(P-P")". 

In the next section, we define properly the BDF functional E'^ and recall its 
main properties proved in [571 UHl HI] • It was noticed in [27] that the minimizer 
of the BDF energy cannot be searched in the trace-class. For this reason, it 
was necessary to extend the definition of the trace in order to give a meaning 
to tr(P — P° ) and to the energy (fT4|) . In this paper, we use all this formalism. 

2.2 The Bogoliubov-Dirac-Fock theory 

We denote by &p{^) the usual Schatten class of compact operators A acting 
on a Hilbert space and such that tr(|A|P) < oo, see, e.g., [50], and by B{9)) 
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the space of bounded operators on Sj. We recall [27] that a Hilbert-Schmidt 
operator A e 62 (^a) is said to be "P^ -trace class if A~^^ = V^AV^ and 
A — = V^AV" both belong to the trace-class 6i{9)a) (but A+- = VlAV" 

and A-+ VlAVl need only be Hilbert-Schmidt). We denote by &^'{^]a) 
this subspace of &2{^a)- We define the T'^-trace of A as 

trpo(A) =tr(A++)+tr(A— ). 

We refer to [37| Section 2.1] for a general definition valid for any reference 
projector and for the useful properties which will be needed in this paper. The 
BDF energy reads p71[^[^ 

£-(Q) :=trpo(I?OQ)-ai?(pQ,i.)-ffi^(pQ,PQ)-f // (15) 

^ ^ J jRe \x — y\ 

where v is the smooth density of charge of a system of extended nuclei, 

and is the free vacuum defined above. We define the BDF energy f on the 
convex set 

Sa |g e ef- (iDA) \Q* = Q, -v°_<Q< rl^ . (le) 

Notice that Q\ is the closed convex hull of the set of operators of the form P — 
V'^ G &2{^a) where P is an orthogonal projector [27l Lemma 2]. Studying the 
BDF energy on Qa will be easier and minimizers will be shown to be extremal 
points, i.e. of the form Q = P — V^. This is a very common technique for 
Hartree-Fock theories [37, . 

Remark 1. Notice that compared to [TU l [TT | HI [7f| . [55] . we have not only 
replaced P° by P° , but also D° by V° in the definition ^ of the BDF energy, 
following the results of p5] . 

Any Q G Qa C &2i^A) has a well-defined integral kernel denoted by Q{x, y), 
such that its Fourier transform Q{p, q) is supported in -6(0, A) x B{0, A). There- 
fore the function Q{x,y) appearing in (|15p is smooth and the charge density 
Pq{x) :— trc4 Q{x,x) is also a well-defined object [57]. In Fourier space, 

p$(fc) = (27r)-3/2 L,/,|<^trc4 {Q{p + k/2,p-k/2)) dp, (17) 

|p-fc/2|<A 

which shows that pg E _L-^(K'^). Introducing the so-called Coulomb space C = 

{/ I D{f,f) < 00}, the linear map Q e ©f" (^Oa) ' — > pq e C n L^{R^,R) is 
v°_ . . 

continuous when &^ {^a) is equipped with the Banach space norm 

IIQIl.pO := ||<3 + +||g^(^^) + llQ ||si($5a) + 11*5^ lls2(flA) + 11*5 ^lle2(flA)' 

as shown in the following useful result, proved in Appendix A. 
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Lemma 1 (Continuity of the map Q i— )■ pq). Assume that < ao < 4/7r and 
A > 0. Then there exists a constant Ca,qo such that 

VO<a<ao, VQeer-(i3A), \\pq\\i^2 + D{pQ, pq)^/^ < Ca^^JQW^.^o . 

(18) 

It has been proved in [27] that is well-defined and bounded-below on Qa, 
independently of A (see also Theorem 1] and [29, Theorem 2.5]): 

VQ e Qa, S^iQ) + ^D{v, v) > 0. (19) 

The proof of (fT^ was itself essentially contained in [3]. If moreover = 0, then 
S'^ is non-negative on Qa, i-e. we recover that is its unique minimizer. 

We shall need to endow Qa with a weak topology for which the unit ball is 
compact. We recall that 6i(i3A) is the dual of the space of compact operators 
[5D1 Thm VI. 26]. It can therefore be endowed with the associated weak-* topol- 

■pO_ 

ogy. This allows to define a weak topology on <Si (^Oa) for which Qn Q 
means Qn ^ Q in 62 (^Oa), 

lim tr{Q++K) = tr{Q++K) and lim tr{Q--K) = tr{Q—K) 

n — >oo n — *oo 

for any compact operator K. It was provecfl in [551 P- 4492] that S'^ is weakly 
lower semi-continuous (wise) for this topology on the convex set Qa, and it 
therefore possesses a global minimizer Q = V- — G Qa where V- is an 
orthogonal projector satisfying the equation 

V- = X(-oo,o] (-Dq) = X(-oo.o] (2?° + a{pQ-iy).l-^- «^^) (20) 

Additionally, if aD{i',iyy^^ is small enough [5Sl Eq. (11)], this minimizer Q is 
unique and the charge of the polarized vacuum vanishes: tr-po (Q) — 0. 



2.3 Existence of atoms and molecules 
2.3.1 Minimization of in charge sectors 

We consider the following variational problem 

E-^{q)= M r{Q) (21) 

Q6Sa(9) 

where Qa('?), the sector of charge —eq, is defined as 

QA{q) Qa, trpo(Q)=(?} 

and q is any real number. 

As recalled before, it is known that the polarized vacuum (i.e. the global 
minimizer of f ) is a solution of E'^{0) when aD{v, vf-^"^ is small enough. But 
in general, it is not obvious at all to prove the existence of a solution to (|2ip . 



^In |28) . the BDF energy is studied on a set 5a slightly different from Qa but the arguments 
used to prove |28l Thm 1] can be adapted by means of Lemma [T] 
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This is because even if the energy functional is wise, the charge sectors Qa(<7) 
are not closed for the weak topology of Qa: a weakly converging sequence might 
loose or gain some charge. We now describe the properties of a minimizer of 
(PT|) if it exists. We recall that m(a) defined in (fTS]) is the threshold of the free 
mean- field operator P". 

Proposition 2 (Self-consistent equation solved by a minimizer). Let be < 

a < 4/7r, v C and g € M. Any minimizer Q, solution of the variational 
problem (PTjl . takes the form Q ~ P — P° — d\(p){tp\, where 

P = Xi-ooA^Q) = X(-oo,H + - ^) * 1/1 • I - (22) 

for some fi G [—m(a),m{a)] and where 

1. if q is an integer, then 5 — 0; 

2. if q is not an integer, then 5 ^ [q] + 1 — q and if is a normalized function 
o/ker {Vq - ^). 

The Fermi level is a Lagrange multiplier associated with the charge con- 
straint and interpreted as a chemical potential. The proof of Proposition [2] is 
left to the reader. It is an adaptation of proofs in [31 [HI IH [571 HH] and of the 
arguments that will be given below for the proof of our other results (see in 
particular Proposition [7]). Notice that when g = A'' is an integer, then ([^ 
means that the last level /x is necessarily totally filled. This is a general fact for 
Hartree-Fock type theories [5] . 

Equation ([5^ is well known in physics. See [JHl EQ- (4)] which is exactly 
equivalent to ([H]) and [371 [H [23 [11 [3D| for related studies. 

Let us assume for simplicity that g = A'^ is an integer. For a minimizer of the 
form ([m and when A^, /x > 0, it is natural to consider the decomposition P ~ 
X(o,m](^q)' where Pvac is the polarized Dirac sea: Pvac := X(-oo , o] (^'q)- 
For not too strong external potentials, the vacuum will be neutral, tr-po (Pvac — 
V^) ~ and therefore X(o,p](2^q) will be a projector of rank N: 

N 

ji=i 

Then Vgipn = e„(/?„, where ei < • • • < ejv are the N first positive eigenvalues of 
Vq counted with their multiplicity. Notice that 

Vq=D + aip^ - ly) * — - a-^-—^ 

1 (Pvac-l/2)(x,y) 
I ■ I \x~y\ 

where p$(a;) :— trc4 (7$(x, x)) — Iv'nl^;)!^- In the first line of ([55)) . the 

Dirac-Fock operator associated with {ipi, ...,ipn) appears, see [H]. This shows 
that the electronic orbitals tpi are solutions of a Dirac-Fock type equation in 
which the mean-field operator Vq is perturbed by the (self-consistent) potential 
of the Dirac sea Pvac — 1/2. Of course, the totally new feature is that these 
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equations have been obtained by a minimization principle (as first proposed in 
[TU] ) while in the Dirac-Fock theory the energy functional is not bounded from 
below. The Dirac-Fock model is thus seen as a non- variational approximation 
of the mean- field model of no-photon QED: the Euler-Lagrange equations are 
similar but the variational structure is very different. In Theorem [3] below, 
we shall prove that the orbitals ((/?!, ...,(^jv) converge to a Hartree-Fock ground 
state [1TJI33] in the non-relativistic limit. 

li N, ^ < Q a, similar decomposition can be applied, 

P = ^'vac - X{p,Q){T^q)- 

When the polarized vacuum is neutral, tr^o (Pvac — Vl) = 0, we obtain 

TV 

P = -Pvac - ^ \(p„){(pn\ 
n=l 

where the minus sign reflects that the orbitals {(pn)^=i describe positrons (up 
to charge conjugation). It holds Vqipn = e-n'fin where e_i > • • • > e_jv are 
the N highest negative eigenvalues of I?q counted with their multiplicity. The 
multiplier fj, is chosen to ensure s-n > > £-n-i- 

Remark 2. To any electronic solution with density v, one can associate a 
positronic solution with density —i' by charge conjugation [271 Remark 8]. 

2.3.2 A dissociation criterion 

The main result of this paper is the following 

Theorem 1 (Binding Conditions & Existence of a Ground State). Let be < 

a < 4/7r, A > 0, V € C and g € M. Then the following assertions are equivalent 

(Hi) for anykeR\ {0}, E^iq) < E^iq - k) + E"{k); 

(H2) each minimizing sequence (Qn)n>i for E'^{q) is precompact in Q\ and 
converges, up to a subsequence, to a minimizer Q of E'^{q). 

If moreover q — N £ 1^ is an integer, then (Hi ) can be replaced by 

{H[) for any K eZ\ {0} , E^iN) < E^iN - K) + E"{K). 

When {H2) holds true, the operator Q is a solution of the self- consistent equation 
(|22p for some Lagrange multiplier fi £ [~m{a),m{a)]. 

Remark 3. Notice that the inequahty E'^{q) < E^{q — k) + E"{k) is true for 
any g e R and any fc G M, as proved later in Proposition [5] 

Remark 4. It will be proved in Lemma [3] below that lim|g|^oo £"^((7) = 00. 
This implies that for any fixed q, there exists a constant M such that |fc| > 
M =^ E^iq) < E^i^q- k) + E°{k). When g = TV > is a positive integer and 
a is small enough (see Corollary H] below for a precise estimate) then it holds 
E^iN) < E'^iN - K)+ E°{K) for all K > N and K < 0. In this case, {H[) 
can be replaced by the more usual condition 

{H[') E^iN) < minlE^iN - K) + E°{K), K ^ 1, ...,N}. 
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Conditions like (Hi) appear classically when analyzing the compactness 
properties of minimizing sequences, for instance by using the concentration- 
compactness principle of P.-L. Lions They are also very classical for 
A^-body Hamiltonians in which the bottom of the essential spectrum has the 
form of the minimum in the r.h.s. of {Hi), as expressed by the HVZ Theo- 
rem [3H im US]. In nonrelativistic Quantum Electrodynamics, such binding 
conditions have also been proved by Griesemer, Lieb and Loss [5Sl[3nj- Notice 
however that, unlike usual HVZ-type results in which a condition similar to {Hi) 
appears only for < fc < iV, here one has to verify that these strict inequalities 
hold for any fc 7^ 0. The reason is that electron-positron pairs can appear. 

For the sake of simplicity, we assume in the following that g = iV is an 
integer. In the next two sections, we provide two regimes in which {H[) is true. 

2.3.3 Existence of a minimizer in the weak coupling regime 

We consider first the weak coupling regime a <C 1 and D := av is fixed (the 
number N of electrons is also fixed) . Our result is the following: 

Theorem 2 (Binding Conditions in the weak coupling regime). Assume that the 
ultraviolet cut-off A is fixed, and that i> G C is such thatk.ei{D^ — i>*\-\~'^) = {0}. 
Then for any integer N G "L, one has 

lim E''/°'{N) = inf trpo {(£)"- i> * | • \-^)Q] . (24) 

Qee^-i^A), -p°<Q<Pl, 

tlpO Q = N 

If we moreover assume that N > and that D Cz C is such that 
(a) a{D^ — D *\ - \^^) contains at least N positive eigenvalues below 1, 
{b) ker(I?o -tD*\ - |-i) = {0} for any t G [0; 1], 

then {H'l) holds in Theorem]^ for a small enough, and therefore there exists a 
minimizer Qa of {N) . It takes the form 

N 

= X(_oo,o] (2?qJ - ^° + X(o,^„] (2?qJ QT + E 10 ('^r I' (25) 

1=1 

2?Q„ = (26) 

where e" < ■ ■ ■ < £% are the N first positive eigenvalues of T>q^ . Finally, for 
any sequence an 0, (v?"", V'at") converges (up to a subsequence) in S)a to 
{(fix, (/37v) which are N first eigenfunctions of — v* \ -\ ~^ and Q™'^ converges 

to X(-oo;0) (D" - j> * I ■ - P- &2{^a)- 

Notice that (5) means that no eigenvalue crosses when t is increased from 
to 1. It is easy to give conditions for which (a) and (5) are satisfied. For instance, 
one can assume that v e Cr\L^{M:'), v > 0, i> ^ Oandthat J^s v < 2/(7r/2-|-2/7r). 
This last constant is related to an inequality of Tix [6OI [61] . 

Remark 5. If is a negative integer, we are able to prove a similar result 
if it is assumed instead of (a) that the spectrum cr(Z?" — D * \ ■ \~^) contains at 
least |A^| negative eigenvalues above —1. 
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2.3.4 Existence of a minimizer in the non-relativistic regime 

Next we consider the non-relativistic regime c 3> 1. For the sake of clarity, we 
reintroduce the speed of light c in the model and we take a = 1 . The free Dirac 
operator is then 

D"{p) ^ ca-p + c^l3. 

The expression of the energy and the definition of the free vacuum V'^ and 
of the free mean-field operator D'^ (which of course then depend on c and the 
ultraviolet cut-off A) are straightforward. To avoid any confusion, we denote by 
E^{N) the minimum energy of the BDF functional. 

In the limit c — *■ oo, we shall obtain the well-known non-relativistic Hartree- 
Fock theory [H] 03], similarly to the non-relativistic limit of the Dirac- Fock 
equations studied by Esteban and Sere in [2^. For a set of orthonormalized 
orbitals (-01, V'Af) S (^^(R^^C^))^, /(V'i,Vi)c2 =<5»j, it reads 

£^Al^):=tr{{-A/2-u.\-\-'h^) + lD{p,^,p,J-]- [[ hfhj^dxdy 

^ ^ J JR6 F ~ V\ 

(27) 

where 7^ — Notice that this model is not posed in i^A but 

rather in the whole space i?^(R"^,C^) since we shall also be able to remove the 
ultraviolet cut-off by taking A — cAq for some fixed Aq > 0. We define the 
Hartree-Fock ground state energy as 

£;Sp(7V):= min £^^p(7v,)- 

Gram^2 ^—Id 

Theorem 3 (Existence of a minimizer in non-relativistic regime). Assume that 
the ultraviolet cut-off is A — cAq for some fixed Aq. Let be v G C D L^{R^, K"*") 
with Jjj3 V — Z , and N a positive integer which is such that Z > A^ — 1. Then, for 
c large enough, {H[) holds in Theorem\^and therefore there exists a minimizer 
Qc for E^{N). It takes the following form: 

N 

Qc = X(-oo,o](2'qJ - ^° + X(o,M.](I^Qc) = Qr + E (28) 

i=l 

and 

hm {E^^iN) - Ngom = E^p{N). (29) 

c — ^00 

Moreover, for any sequence Cn — > 00, (c/?J", c/?^) converges (up to a subse- 
quence) in iJi(M-\C'^)^ towards ip = (j^^ iP £ H'^{R^ ,C'^)^ , and where 7^ is 
a global minimizer of the Hartree-Fock energy (|27p . 

It is proved in Lemma [T9l that go(0) = min(T(|I?°|) is the threshold of the 
self-consistent free Dirac operator I?" for c large enough. 

The rest of the paper is devoted to the proof of Theorems [1] [5] and [31 

3 Preliminaries 

3.1 Behavior of E^iq) for \q\ > 1 

We give conditions which prevent the appearance of electron-positron pairs in 
minimizing sequences. 
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Lemma 3. Assume that < a < 4/7r, A > and u ^ C. Then one has 

(1 - aTi/A)m{a)\q\ - -D{v,v) < E'' (q) < .go(0)|g| (30) 

where (?o(0) > 1 is defined in (jlOp and 771(0;) is defined in (jl3p . Hence it holds 
\\m.\q\^oo E'^ {(l) =^ 00. For any fixed q € M., there exists an M depending on q, 
a, V such that \k\ > M =^ E^lq) < E^l^q - k) + E°{k). 

Proof. For the right hand side of (jSn)) . let us fix some orthonormal system 
{ipi, ■■■,ip[q]+i) of smooth C^-valued functions with compact support in the 
Fourier domain. We introduce the following 

^^. = Span{7'>^ 7 = l,...,[g] + l}, = ) = 

Note that for A small enough W\ is a subspace of Sj'^ of dimension [q] + 1 since 



^V9o{\p\r + 9i{\p\y 



Let us choose an orthonormal basis ((^^, c^j^j^-^) of Wx. The r.h.s. of ([50]) is 
then obtained by taking a trial state of the form 

/ [<?] 

where e = 1 if (7 > 0, and e — —I otherwise, and by taking the limit A — > 0. To 
prove the lower bound in (|30p . one uses that 4J for any Q E QaI?), 

S^UFiQ) > (1 - an/A) trpo (I?°Q) - ^Dii,, ly), 

trpo (pOQ) = tr(pO|(Q++ _ Q— )) > m{a) tr(g++ - Q— ) > 77i(a)|g|. 

□ □ 

Corollary 4. Let be < a < 4/7r, 7/ G C and A > 0. Assume that N is a 
non-negative integer and that 

(go(0) - 777(a))iV + a (^m{a){N + 2) J + R^!^^ < 2777(a). (31) 

Then E''{N) < E''{N - K) + E^{K) for any integer K satisfying K > N or 
K < 0. Therefore, in this case the HVZ-type condition {H[) in Theorem]^ can 
be replaced by the more usual one 

{H'{) E" (N) KminiE" {N ~ K) + E"{K), K^1,...,N}. 

The proof of Corollary |4] is left to the reader. When m{a) = 5o(0) (which is 
true when a ^ 1, see Lemma fTO)) . (|5T1) can be replaced by the stronger condition 
a((7V + 2)7r/2 + i:)(7^,7/)) < 4. 
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■po 

3.2 Approximation by finite-rank operators in 

Proposition 5 (Approximation by finite-rank operators). The set consisting 
of the operators Q which satisfy 

1. Qe QAiq); 

2. Q = P — + 7 where P is an orthogonal projector and j is a finite rank 
operator such that < 7 < 1, P7 = 7P = 0; 

3. Q has a finite rank; 

is a dense subset of Qa{(i) for the strong topology of &i . 

Proof. The proof relies on a useful parametrization of the variational set Qa, 
presented and proved in Appendix B, Theorem[6l This result itself is a general- 
ization of a reduction in the case where Q € Q\ is a difference of two orthogonal 
projectors, see Theorem [SJ By Theorem [51 any Q € Qa{N) can be written 
Q = Ud{V1 + j)U-D - V° where 7 € &i{^a), [j,V1] = 0, D e e2(iOA) 
and Ud = exp{D — D*). Moreover 7 = 7+ — 7^ where < 7"'' < V^. and 
< 7~ < V?_. Clearly we can find sequences {in} of finite rank oper- 

ators such that Dn D in &2{SiA), In ~^ 7^ ™ &i{J^%S)a) as n ^ 00, and 
additionally tr(7+-7-) = tr(7) for any n. Then Ud^VIU-d^-UdVIU-d ^ 
in 62 (-$3 a)- We know from [57| Lemma 2] that 

is an integer. Thus tr^o {Ud„P-U-d„ —V^) is constant for n large enough and 

is a sequence of finite rank operators which converges to Q in QaI?)- Assume 
now that 

it = J2^t\vt){^t\ and 7,7=E^^"l^r>(^n 

i=l i=l 

with < Xf <1 and e P^-^a- We then introduce 

^^=^--El^r>(^n+ E \v^t){^t\ and y,^:=vl+jt~ln~Pn- 



Then Q„ — Ud„ (Pn + l'n)U-D„ — P- = Pn + 7ri — P- satisfies the assumptions 
of the Proposition. □ □ □ 

Corollary 6. There exists a minimizing sequence {Qn)n>i of E'^{q), satisfying 
the three conditions of Proposition [5l 

As usual in Hartree-Fock type theories [37l [3l [6] , we now prove that mini- 
mizing S'^ in the convex set of states in Qa having charge —eq is equivalent to 
minimizing on extremal points only. 
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Proposition 7 (Licb's variational principle). Let beO < a < 4/7r, A > 0, G C 
and q G M. One has 

E''iq)^mi{riQ) \Q = P-V°+{q-[q]m{^\eQA, 

P^ = P = P*, trpo(P-7'0) = [g], 7/;Gi3A, = O}. (32) 

Proof. We use a well-known method for the study of Hartree-Fock type theories 
[371 [31 [B]. The inequality < in ([3^ being trivial, it suffices to prove the converse. 
To this end, we consider by Proposition [S] a state Q — P — V^ + 7 G Qa(<z) such 
that = P = P*, 7 is a finite-rank operator with < 7 < 1, P7 — 7P — 0, 
trpo Q ^ q, and £'^{Q) < E'^{q) + e. If 7 has rank greater than 1, then one can 
find two orthogonal eigenfunctions xi ^-nd X2 of 7 corresponding to eigenvalues 
71 and 72 in (0; 1). Next one can compute the energy of Q+t{\x2){X2\ — \xi){Xi\) 
which belongs to QA{q) as soon as 72+^ G [0; 1] and 71 G [0; 1]. One sees that, 
depending on the sign of {T>qX2, X2) — {'DqXIi Xi), we can decrease the energy 
by increasing or decreasing t until 72 -I- i G {0; 1} or 71 — t G {0; 1}. Doing so, we 
obtain a new state Q ^ P — P£ + 7 where P is an orthogonal projector and 7 
is such that < 7 < 1, P7 — 7P = and rank(7) < rank(7) — 1. Iterating this 
process, we can eliminate in finitely many steps all the eigenvalues in (0; 1) of 
the finite-rank operator 7, except possibly one, and decrease the energy at each 
step. We end up with a state Q' = P' -V°+ A|'0)(V'| where A G [0; 1), PV = 
and £'^{Q') < E'^{q) + e. It then suffices to use jUl Lemma 2] which tells us 
that ti-pa (P' — Pi) is always an integer, to conclude that X = q — [q]. □ □ 

3.3 HVZ-type inequalities 

We use the density of finite-rank operators to prove HVZ-type inequalities. 

Proposition 8 (HVZ-type inequalities). Let be < a < A/tt, A > 0, ly E C 
and q G M. Then one has 

E^iq) < min {E" {q ~ k) + E°{k), fc G M} . (33) 

// moreover E^{q) — E^{q~ k) + E''\k) for some fc 7^ 0, then there exists a 
minimizing sequence {Qn)n>i of E'^{q) which satisfies Qn ^ Q in Q\ with 
tr-po (Q) = q ~ k, and which is therefore not precompact in Qa('?)- 

Proof. Notice that E'^{0) = by [571 Theorem 1], so there is nothing to prove 
when k = 0. Let us fix some k Q and consider two states Q, Q' G Qa such 
that £°{Q) < E°{k) + e and S^iQ') < E''{q - fc) e, e > 0. Using Proposition 
[5l Proposition [7] and Theorem [3 we can choose Q of the following form: 

Ni N2 K .2 

Q = A|^)(V'|-f E l5™)(5m| + ET:rA^(l^^)(«^l-l"^)("'l) 

n— 1 m— 1 i—1 * 

{\ui){v,\-rW){u,\). (34) 



where A = /c - [/c] G [0; 1), , A^i - TVs = \k\. In ([Mj) {fi)^=x U (i;,)^i forms an 
orthonornial set of P^^iOA and {gi)^Ji U (ui)|£i is an orthonormal set of P^^da- 
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We may choose the operator Q' of a shnilar form (with fmictions V'', f'm g'm' ""ii 
vl and real numbers A', A-). 

The free minimization problem E^(k) being translation- invariant, we have 
S^inQ'Tf) — £^{Q') for any i G M, where Tt is the unitary operator Ttf{x) = 
f{x — te), e being a fixed unitary vector in M.^. Notice nQ'r^ € Qa, since both 
and are translation- invariant. Since limj^oo {Ttf,g) = for any f,g G 
Sja, we can find for f 3> 1 by the Gram-Schmidt orthonormalization procedure 
two orthonormal systems {fn)n=i ^ i^Dj^i in the orthogonal of span{/^, v'^} 
inside 'P+i^A, and {gD^Ji U in the orthogonal of span{g^„, d-} inside 

V-^A which are such that lim^^oo ||/^ — rt/„ jj^^^ — and a similar properties 
for all the other functions. Substituting in (|33|) . this defines us an operator 
Q* £ Qa such that hmj^oo ||Q* — TtQrj* Ig^^^^-j = 0. Moreover, we have by 
construction + Q' e Qa and tr^o (Q* -I- Q') — q. It can be seen that 

E-'iq) < E^Q' + Q') = E^Q') + £"(Q) + ot^^{\) 

<E''iq-k)+E°{k) + 2e + Ot^^.{l), (35) 

which, passing to the limit as i ^ oo and e ^ 0+, ends the proof of this first 
part. If E'^{(p) — E'^^^{lp) + E''(0) for some A: ^ 0, then one constructs a 
non-compact minimizing sequence by the same argument. □ □ 

Corollary 9. Let be < a < A/n, v eC and A > 0. Then the map q ^ E"{q) 
is uniformly Lipschitz on R; 

\E^{q)-E''{q')\<g^{Q)\q-q'\. 

Proof. This is an obvious consequence of and ([30]). □ □ 

4 Proof of Theorem [1] 

We prove that {Hi) implies {H2), as the converse was shown in Proposition [H 
Step 1 : Reduction to the HVZ condition for integers when q = N G "L. 

Lemma 10. Let be Q < a < 4/7r, v G C and G Z. Then for any K eIj, the 
function k 1— s- E'^{N — fc) + E'^{k) is concave on [K, K + 1]. Therefore, if {H[) 
holds, then so does (Hi). 

Proof. We prove that the function k 1-^ E'^{k) is concave on [K, A'+ 1], the rest 
being obvious. This means that we prove for any /ci, fc2, /i G [0, 1] 

E''{K + nki -f (1 - ^)fc2) > fiE-'iK + fci) + {1- n)E''{K + ka). (36) 

To this end, let us consider like in the proof of Proposition [HI one state Q 
which satisfies f (Q) < E''{K + fiki + (1 - fJ.)k2) + e for some fixed e > 0. By 
Proposition [3 we may take Q — {fiki + (1 — fi)k2)\ip){^j\ + P — where P is 
an orthogonal projector satisfying Pip = and tr^o (P — V'^) = K. We have 

E-'iQ) = {fiki + (1 - /i)fe)(2?[P-pO]V', ^) -f E'iP - Vl) 

= iiE^P -Vl + fci|^)(^|) + (1 - fijE-^iP -Vl + 
> fiE-'iK + ki) + (1 - ^l)E''{K + fcz) 
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where we have used that an electron never sees its own field. This ends the 
proof, taking e — > 0. □ □ 

Step 2: A necessary and sufficient condition for compactness. We now assume 
that g S M and prove that conservation of the charge implies the compactness 
of minimizing sequences. 

Lemma 11. Let be < a < 4/7r, A > 0, G C and q £ M.. Assume that 
{Qn)n>i is a minimizing sequence of E'^(q) such that Qn Q for the weak 
topology of Qa. Then Qn Q for the strong topology of Qa if and only if 
tr-po {Q) = q. 

Proof. By assumption, we have lim„^oo {Qn) = E'^il)- If tr-po (Q) — q, then 
Q is a minimizer for E'^{q) since is weakly lower semi-continuous. Also we 
have [28] hm„^oo \\pq„ - PqWq = and lim„^oo F{Qn) = F{Q) where 

F{Q) := tr^o (V^Q) ff ^^^dx dy. 

2 77r6 \x - y\ 

Let us first prove that Q„ — > Q in S2(^a), or equivalently tr{Qn) — > tr(Q^). 
To this end, we argue like in the proof of (23 Theorem 1] and consider two 
smooth functions x, ^ G C°°([0; oo), [0; 1]) such that + = 1 and 

{= 1 when X G [0; 1] 

e [0;1] whena;e[l;2] 
= when x > 2. 

We then define xr{^) •= x{\^\/R) and £,r{x) :— ^(|a;|/i?) for x G R^. As shown 
in [3S1 page 4494], one has 

F{Qn) = F{Q) + tyi\V°\UQt^ - 

_arr U-?\Qni^:y)\\ ^,y^^H (37) 

2JJk6 \x-y\ 

where satisfies limfl^^oo linisup„_o2 |e^| =0 (see the details in [21] )■ Local- 
izing the inequahty Q^ < — Q~~ and applying Kato's inequality < 
7r/2|V| like in [3S] together with |V| < \D"\ < \V°\ as shown in [H], we get 

^ J J«6 \x-y\ 
> (1 - a7r/4) ir{\V°\iRQliR) > (1 - a7r/4) tr{^BQl^R)- (38) 

By ((311), ((H) and i^(Q„) F{Q), this proves that 

lim limsuptr(^ijQ^^ij) = 0, 

R 'OO fi — 

when < a < 4/7r. On the other hand 

tr(g2 ) = trixRQlxR) + tri^RQliR) 

and 

lim tr{xRQ'ixR) = ^'^{xrQ'^Xr) 
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for any fixed i?, by the local compactness of Qn. This implies 

lim tr(Q2 ) ^ tr(Q2) 

n — >oo 

or equivalently Qn Q in &2{^a)- 

It remains to prove that Q^^ and Q^^ converge strongly in the trace-class. 
We use the continuity of the exchange term for the Hilbert-Schmidt norm 



' ^ dxdy<-tY{\W\R^) <—\\R 



|2 



\x-y\ - 2 ^^^1^1"^ ' - 2 "''"©2(JjA)' 
due to Kato's inequality and the cut-off in Fourier space, to infer 

lim // \9^d.dy^[[ ^^dxdy. 
"-'°°7Jk6 \x-y\ JJr6 \x-y\ 

This proves 

lim tr(|2?0|(Q++ - Q") = tr(|2?°|(Q++ - Q") 

n — *oo 

and therefore (5++ and Q in &i{Sj\). □ □ 

We finish the proof of Theorem [T] by a contradiction argument: we assume 
that (Hi) holds and that there exists a minimizing sequence {Qn)n>i of E^^q) 

which is not precompact for the strong topology of (-^a)- Since £^(Qn) 

converges to E^{q), (Qn) is a bounded sequence in 6^ {^a) and we may there- 
fore assume, up to a subsequence, that Qn Q for the weak topology of 

&i (-Qa)? and that Q„ Q. By Lemma this is equivalent to assuming 
that trpQ (Q) ^ q. We now write tr^o (Q) = q — k with k E M. \ {0} and prove 

that this would imply E^iq) >E''{q-k) + E'^(k), which contradicts {Hi). 

As usual in HVZ or concentration-compactness type arguments, the rest of 
the proof now proceeds by decomposing the sequence (Qn) into a compact part 
converging strongly to Q and a non-compact part escaping to infinity with the 
charge k. The localization of Q is complicated by the constraint appearing in 
the definition of the variational space Qa((7), and the fact that our states are 
not bounded in the trace-class. 



Step 3: The localization operators. We introduce 

Xr VIxrVI + VIxrK 

and Yr which is the unique non-negative operator satisfying X]^ -f = 1^5^. 
Recall that the function XR has been defined in the previous step. A crucial 
fact will be that both Xr and Yr commute with P". The next three lemmas 
summarize useful properties of Xr and Yr. 

Lemma 12 (Continuity of the localization maps). There exists a constant C 
independent on R and A such that, for any Q E &i (-$3a)j 

\\XrQXr\\^..po_ + \Yr,QYr\\i,v«_ < C IIQIi.po . (39) 

If moreover Q g Q\, then XrQXr and YrQYr also belong to Q\. 
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Proof. We notice that {Xj^) and (Y/j) are uniformly bounded in Soo(-^a)- In- 
deed, using < xr 1^ 1; one sees that \\Xfi\\ < 1 and ||Yff|| < 1. Therefore 

W^RQXHWe,^^,^ < WXnf IIQIIe.(^.) < IIQIe.c^.) > 
and, using that commutes with , 

\\iXRQXR)—\\^^^^^^ = WXhQ—XrW^^^^^^ < ||Q— 

by the same argument as above. Using the same idea for YrQYr, one obtains 
(jS^ . Let us now prove that if Q satisfies the constraint —V^ < Q < 
then so does XnQXj^, the argument being the same for YrQYr. We have 
XrQXr < XrVIXr = [VlxRVlf < VIxrVI < Vl- A similar argument for 
XrQXr > —V^ allows to end the proof of Lemma [T2l □ □ 

Lemma 13 (Limits as i? oo). One has 

lim \\Xr-xr\\=0, lim ||yj^-ei?||-0. (40) 

R — >oo R — ^oo 

■p" 

Moreover, if Q (z &i {^a), then 

lim WXrQXr-QI.^o = lim \\YRQYRi,.^„ =0. (41) 

Proof To prove ^ for Xr, we compute Xr-xr^ [P^,Xr]'P°^ + [P+, Xr]'P+ 
which tends to as i? ^ cx) by ^5] Lemma 1] and the proof of Theorem 3 in 
pS] , This clearly implies that limi^^oo H-'^l; ~ X?j|| = 0, since \\Xr\\ < 1 and 
WxrW < 1- We now use the fact that the square root is operator monotone to 
deduce [7, Thm X.1.1] 

\\YR-u\<m-eRf" = m-x'Rf' 

which proves (|40|) for Yr. By the uniform boundedness of {Xr)r in 6oo(-J3a), we 

can prove ((41|) for Q in a dense subset of &i (S^a), like finite-rank operators. By 
linearity, it suffices to prove (jUJ for a state of the form Q = Using now 

(|40|) . it remains to prove that XrW WXR — converges to in 6i(^a), 

which is a triviality since X-rV' — > ?A as i? — > oo in Sja, by Lebesgue's dominated 
convergence theorem. We argue similarly for YrQYr. □ □ 

Lemma 14 (Compactness for a fixed R). For any fixed R, Xr and 1 — Yr belong 

to Si(i3A) cind are therefore compact. The map Q i— s- XrQXr is also compact: 

■po 

if Qn Q for the weak topology of &-^^ (S^a) , then XRQnXR — ^ XrQXr 
strongly in 6i {S)a)- The same holds if Xr is replaced by 1 — Yr. 

Proof. We use the Kato-Seiler-Simon inequality (see and ^Tj Theorem 4.1]) 
\\fH^)9i^)\\e.<C\\fh2 bh^ 

to obtain 

Il^i..|l6,(^.) < ll^-V^ire.(«.) + ll^+V^llli^.) ^ 2C|i?(0, A)|i? f x 
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which is finite for each fixed R. Notice that < 1 — ^7? < ^^^r = ^r ^ ^R ^^'^ 
therefore 1 — Yr is also in &i{Sja). Let us now prove that when Q„ ^ Q for the 

weak topology of 6^ (^a), then XrQuXr -> XrQXr strongly in 6^ (i^A)- 
Since Xr is in Si(^a) and by density of finite-rank operators in the trace-class, 
it sufhces to prove that \(p){(p\ Qn\'4'){'^\ l¥')(</j| Q for some ip and 

ij] in i3A- This is a consequence of Qn — ^''^{Qn\ip){v\)\v){''P\ ^^nd 

lim„_»oo tr((3„|V)('Pl) = tr(QIV')(</'l) by the weak convergence of Qn in S2(i3A)- 
The proof is the same for 1 — Yr. □ □ 

Step 4: Conclusion. Let us recall that we argue by contradiction: we as- 
sume that {Hi) holds and that [Qn] C Qa is a minimizing sequence for E'^{q) 
satisfying Qn — ^ Q with tr-po [Q) = q — fc, fc 7^ 0. We fix some R and compute 

trpo (I?0Q„) = trpo (I?"(Xi^g„Xfi)) + trpo {V"{YrQ„Yr)) 

+ tr([Xfl, |i?°|](Q++ - Q,T")^fl) + tr([rfl, |i?°|](Q++ - 0--)r^), 

where we have used that Xr and Yr commute with Using now 
I triiXR, |I?"|](Q++ - Q--)Xr)\ < C \\[Xr, 

due to the fact that {Xr) is uniformly bounded in B{S)\) and {Q^^ — Qn ) is 
uniformly bounded in 6i(^a), we obtain for some uniform constant C 

trpo {V'>Qn) > trpo {V°Q) + tipo (i?°(ri?g„rfl,)) 

+ trpo {V^{XR{Qn - Q)Xr)) + tr-po {V"{XrQXr - Q)) 
-C{\\[Xr,\V''\]\\ + \\[Yr,\V'>\]\\). (42) 

To treat the direct and exchange terms, we shall need the following 

Lemma 15. Assume that {Rn) is a sequence in &i (-^a) with Rn ^ as 
n — > 00, for the weak topology of (^a)- Then, for any fixed R, 

lim \\PR^ - PYrR^yAc = 0, (43) 

n — >-oo 

hm // '(^"'^f^"^f)(-'^)'^..d, = 0. (44) 

Proof. Recall that the exchange term is continuous for the 6^ {^a) topology. 
We write 

Rn - YRRnYR = -(1 - YR)Rn{l - Yr) + i?„(l - Yr) + (1 - YR)Rn. 

By Lemma [Ml (1 — YR)i?„(l — Yr) converges strongly to in the 6^ {^a) 
norm. It therefore suffices to prove ([33]) and (HH) with i?„ — YRRnYR replaced 
by — Yr) and (1 — YR)Rn. The operator (1 — Yr) being trace-class by 

Lemma [HI we can prove (|43l) and (|44l) with i?„ — YRRnYR replaced by S'„ := 
Rn \ip){ip\ e 6i(i3A) for some G i^A H ^^(M^) C L2(r3^(;-4)^ j^^^^g 

Ps„{x) ^ / trc4 (i?„(a;,y)(^(y)v3(a;)*)d?/ 
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and therefore 

||ps„|lii < // \Rn{x,y)\\ipix)\\ip{y)\dxdy. 



Thanks to the cut-ofF in Fourier space, we can assume, up to a subsequence, 
that Rn{x, y) converges uniformly to on compact subsets of M^. Since we have 
assumed that (p S L^{M.^), we conclude by Lebesgue's dominated convergence 
theorem that lim„_,oo ||ps„|lii = 0, which itself implies lim„_,oo lAs^iioo = 
and lim„_,oo ||/Os„|lc ~ thanks to the cut-off in Fourier space. 
We use the same argument for the exchange term: 

\Snix,y)\^ ^^^^^ ^ ff 1^ .m2|.„/ m2 A, ,2 ... 1 



F - y\ 



dxdy< // \Rn{x,z)\'\ipiz)\^[\ip\^^--]{x)dxdz 



which converges to as n — > oo by the local compactness of (i?„). This ends 
the proof of Lemma [T5l □ □ 

We are now able to finish the proof of Theorem [1] We write 

D{pq,,,PqJ = D{pq,Pq) + D{pY„(Q„-Q)Ya, PYa{Q„-Q)Yn) + ef H 

> D{pq,Pq) + D{pYj,Q„Yn,PYnQ„Yn) + ef (") - Ci \\pYnQYn\\l ' (45) 



ff \QA^,^,y 

J Jjj6 \x — y\ 

JJffie \x-y\ JJm.6 \x-y\ 

^ ff \Qix,y)\' , , _^ ff \{yB.QnYR){x,y)W 

< // — —dxdy+ // . . dxdy 

JjRe \x-y\ JJs,6 \x~y\ 

^ R( ff \YRQYRix,y)W , .... 

+ £2 n) + C2 . . dx dy, (46) 

JJr6 \x-y\ 

and 

DipQ,^,iy)^D{pQ,iy) + 63(71). (47) 

In (|T7| and Ci and C2 are uniform constants (we have used that 

■po 

(Qn) is bounded in " (^a)), and £3(71) := D{pq^^q,v), 

ef(") IIpq„-qIIc- \PYn{Qr.-Q)Yj\c ^ '2D{pQ^^Q, pq), 



\{Qn-Q){x,y){' 



\x - y\ 



■dx dy 



\YniQn-Q)YR{x,y){' 



\x - y\ 



■dx dy 
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trc4(Q» - Q){x,y)Q{.x,y)] 

\x - y\ 



dx dy (48) 
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satisfy, for each fixed R, lim„^oo ("■) — lini„_>oo ^l^l"-) = lini„^oo ^ain) — 0, 
by Lemma [T5l applied to i?„ = Qn-Q- By gH), gS]), (gT]) and we obtain 

f "^(Qn) > S^iQ) + £"{YRQnYR) + ef (n) + - (n) 

+ trpo (I?0(X,j(Q„ - g)Xfl.)) + trpo (I?0(Xi^QXi^. - Q)) 



C \\[Yn, \V°\]\\ - C \\[Xr, \V°\]\\ - Ci IpYnQYn 



~CJ[ M^dxd,. (49) 

Let us remark that q ~ ii-pa {Qn) — tr^o {XjiQnXB) + trpo {YjiQnYR) where 
we have used [Xr.V'^] = 0. In ([^5]) . we can estimate 

£'^(g)+£"(r^Q„lfl) > E%q-k)+E\iYT,<^{YRQ,,YR)) 

= S''(g-fc)+S"(g-trpo(Xj^Q„XK)). 

Passing to the hmit n ^ oo in (|49p with i? fixed and using Lemma [TJ] together 
with the continuity oi E'^il) a-s stated in Corollary [3 we get 

E'iq) >E''{q-k) + E°{q - tr^o (XrQXr)) + trpo [V^iXRQXR - Q)) 

-C\\[Xr,\V''\]\\-C\\[Yr,\V''\]\\ 

r II i|2 r ff ^WflC^i^lP. . ,,nA 
- iPi-RQyR ic ~ / / . . dxdy. (50) 

J Jrb F — 2/1 

Note that 

hm = hm ||[XK,|I?"|]|| = hm || [Ffl, || 

= hm =0 

by Lemma [13] together with [23 Lemma 1], [3^1 Proof of Thm 3] and the bound- 
edness of 2?° on ^a- Passing to the limit as i? ^ oo we eventually obtain from 
Lemma [T3I 

E-'iq) > E^iq -k) + E\q - tr^o (Q)) = E'' {q - k) + E^k). 
This contradicts {Hi) and ends the proof of Theorem [T] □ 

5 Proof of Theorem [2] 

Note that when v & C, the essential spectrum of — P * | • is the same as 
the one of . Assuming that ker(Z?" — j/ * | • |~^) — {0}, we denote by (A^)i>i 
the non-decreasing sequence of eigenvalues of D° — P * | • in (0, 1), counted 
with their multiplicity. In case there is a finite number ig of eigenvalues or no 
eigenvalue at all (zq = 0) in (0, 1), we then let A^_|_;^ = 1 and A,^ = 1 for i > io + l. 
We use the same type of notation (A~)i>i for the non- increasing sequence of 
eigenvalues in (—1,0), with AT7^-^ = —1 in case there is a finite number of 
(possibly no) eigenvalues in (—1,0). We notice that ker(_D'^ — v |^^) ~ {0} 
implies that there exists some constant k > such that 

|i?0-//*|.|-i|>K. (51) 
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step 1: Study of the linear model. We start by computing the value of the 
infimum in the right hand side of in terms of the eigenvalues (A^) of 

Z)° — D * I • 1^^, and the charge of the vacuum. We define 

r(iV) ^ inf trpo {(£)"- I ■ r^)*?} • (52) 

Qeef-(i5A): -p°<Q<p!1, 

trpo Q=N 

Lemma 16. Assume that v ^ C is such that k.eT{D^~iy*\ -l^^) = {0}, and denote 
the charge of the non-interacting (Furry) Dirac sea by qo := trpo (P — P^) 

where P := X(-oo:0)(^° — v *\- \^^). Then one has for any iV G Z 

\N~qo\ 

r(7V)^trpo z>H-r')(^ E l^'^l (53) 

where e = sgn(Af — qo)- 

Proof. After a change of variable Q ^ Q — {P ~ P^), and by [27l Lemma 1], 

r (TV) - trpo {(D" - * I ■ r^) - -P°)} 

inf trp{(i50-P*|-r')Q}- (54) 

Qesf (s5a), -p<q<p, 

tip Q=N-qo 

By a simplified version of Proposition [71 one sees that the infimum of the r.h.s. 
of ([54]) can be restricted to states Q which are a difference of two projectors: 
Q = P — P (recall qo E Zhy [27l Lemma 2]). By Theorem [5] proved in Appendix 
B with n = P, there exists two orthonormal basis igi)fji U (ui)i>i and {fi)^J,i U 
{vi)i>i respectively of PS)k and (1 — P)S)a, and {k) G £^(M) such that: 

Ni N2 j2 

n— 1 m— 1 2>1 ^ 

+ Et^(I"»)(^^i + i^^)(^^i)' (55) 

with N — qo = iVi — A^2- The following inequality gives the lower bound in ([55]) : 

trp{(D"-z>H-r')Q}>E^«^-Ev 

The proof of the upper bound is left to the reader. □ □ 

Step 2: Upper bound. To obtain an upper bound for (|24p . we first fix some 
> and a state P — P^ + 7, 7 being a finite-rank projector commuting with 
— * I • |~^, which satisfies 

tr(7) =N-qo (56) 

and 

7 > and 7P = P7 = if iV - go > 0, 
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7 < and 7-P = ^'7 = 7 if ^ - 90 < 0. (57) 
We additionally assume that 

ivpo {{D° |-i)(P - P" + 7)} < F{N) + 77. 

If for the considered charge N , then only eigenvalues appear in formula (|53p. 
one can of course choose 77 = and 7 to be the projector on the space spanned 
by any chosen eigenvectors associated with the (A-). However, if Af = ±1 for 
some i, then a minimizer does not necessarily exist and we can only take an 
approximate one as expressed above. Then, we recall that [IHl [35] 



0{a) 



75U 



0{a) 



(58) 



This in particular implies that the spectrum oiTf' — v ■ 



of £)0 



as a ^ 0. Therefore, ker(I?" - 
|pO__ - ^ I . |-i| > 



^ converges to the 
1^1) = {0} and 

(59) 

for a small enough. We shall consider a trial state of the form X(-oo;0) (2^" ~ 
* I • l"""^) — + 7q, where 7^ is a projector converging to and having the 
same rank as 7. We assume moreover that it satisfies the same properties (j56p . 
© and jlTl) as 7, with P replaced by Pa := X(-oo;0)(2^° - P * | • 1"^). This is 
possible, since P^ ~ P in B{f)fO as a 0. We then introduce 

Qa := Po.-V'L^ X(-oo;0)(2?° - ^ * I • r') - ^- (60) 

(recall that I?'^ and P'^ implicitly depend on a) and prove the following 

Lemma 17. The operator Qa defined in (|60p satisfies: 



lim 

a— >0 



Q„-(p-p°; 



62(flA) 



= 0. 



Moreover, one has for a small enough 

trpo(Q„) =trpo(P-P") = 90. 
Proof. Consider first the simplified case where (f := —D> 



(61) 
(62) 

-1 e L2(M3). Write 



(P-P^) = 



1 

2^ 



dr] 



1 



+ir]^V" - P * \ ■ + it] 



1 

2^ 



+00 



1 



1 

2^ 



4-CXD 



+00 



1 



177 



£)" +iri^V° -v*\-\-^ +iri 



X)0 



IT] 



1 



1 



+ 17] - P * I • |-i +177' 



(63) 




V(1 + 772)(k2/4 + 772) 



(64) 
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-if 



< C 



by the Kato-Seiler-Simon inequality |561I57| and a similar estimate for the second 
term of to obtain 



<C\\V" -D°\ 



B{SOA 



= 0(a). 



To treat the general case, it then suffices to approximate D appearing in Qa 
and P — P° by a function such that * | • G L^, uniformly with respect 
to a. This can be done by using the method of [SB] 123 • '^^^ be shown that 
there exists positive constants Ci and C'2 independent of a (but which depend 
on D) such that, for any p satisfying — < Ci, 



(X(- 



oo;0) 



62(flA) 



<C2\\i>-iy\ 



C ' 



<C2||i>-i/|| 



c • 



62(^1/1 

Taking for instance Pe{k) — S'(A;)ln;|>e this allows to end the proof of ([61]) . 

To end the proof of Lemma [TTl one notices that since Qa is a difference of 
two projectors, 

Ql = VlQo.Vl - VlQe.V'L (65) 

and similarly (P - pO)^ = pO(P - pO)pn - P^{P - P^)P"_. Therefore, (EH) 
implies that 



lim 

a->0 



P'^QaP^+P°OaP° 



Si(i5A) 

= \\Pl{P - P°)Pl + P°{P - P°)P'1\ 



from which we infer lima^o tr-po (Q^ ) — trpo (P — P^ ) . Then 
since both are integers by \2I, Lemma 2]. □ 



II6i(S3a) 

is proved 
□ 



For a small enough, we deduce from (f62|) that tr-po (Qa +^a) — N and that 

Qa + 7a is an admissible trial state. Note that 7^ is bounded in hence 

lim^^o ^BDF(7a) ~ tr{(P)° — u * \ ■ |~^)7}, as the direct and exchange terms 
vanish in the limit (they are multiplied by a). 

Since (Qa) satisfies (|65l). it is bounded in Qa- This imphcs that (pg^) is 
uniformly bounded in C by (|18p , hence 



V I OL 



BDF(Qa + 7a) = ^'^V'l [T^^Qo) - P>{v , PqJ + ^BliFila) + 0(a). 



Using one more time the fact that Qa is the difference of two projectors, one 
deduces that tr^o CD^Qa) = ir{\V^\Ql) which converges to tr(|D°|(P-P°)2) = 

trpo(D°(P - P")) by Lemma [TTl and (f55| . Now, using dHH) and the fact that 
{Pq ) is uniformly bounded in C, we obtain lima^o ^'(t^j Pq ) = D{i>^ pp_pa). 
By [271 Lemma 5] which ensures 



trpo (i?"(P - P'_)) - D{9, pp_po ) - trpo {{D° - 9 ^ \ ■ \-') (P - P^)} , 
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we have proved that 

\im r^i,%{Q^ + 7.) - trpo {{D'' - D * \ ■ \~^) {P - P"_ + 7)} , 

a — *U 

which means hmsup^_^Q i?g(j^(A^) < I'^{N) + rj for any ?7 > 0. 

Step 3: Lower bound. To prove the lower bound, we consider for any fixed a 
a state {Qa) satisfying 

4df(Q«) <^'^"W + «- trpa(Q„) = iV. (66) 

By Proposition [3 we may moreover assume that Qa — Pa — "P- where Pa is 
an orthogonal projector. Let us show that the sequence {Qa) is bounded in 
'Sii^A)- To this end, we first notice that E'^^°'{N) is bounded from above by 
the previous step and therefore, by and Kato's inequality for the exchange 
term, we obtain the bound 

(1 - air/i) tr(|I?0|Q2 ) <c+ ^D{D, P), 

which proves that {\/aQa) is bounded in &2{^a)- By Kato's inequality, this 
means that the exchange term satisfies 



a 



11 \^^^d.dy < ^tr(|pO|Qi) < + 

JJvfi \x-y\ 4 



2yyRe \x-y\ ^' 4(l-a^/4) 

and it is thus uniformly bounded in a. Therefore, one has 

trpo (P°Q„) - D[V, pQ^ ) + ^D{pQ^ ,PqJ<C (67) 

for some other constant C independent of a. By 27, Lemma 5] we have 

tTMV°Qa)-DiD,pQj = trpo{(I?o„j^^|.|-i)Q^} (68) 

= tYpJiV" -D*\-\-^){Qa-Qa)} 

+ trpo{(I?"-P*|-r')Qa} (69) 
> tr{|I?" - p * I . |-i|(Q„ - g„)2} _ a (70) 

In dMl), we have inserted Qa = Xi-oo-fi)^^" - v * \ ■ \-^) - = Pa - Vl used 
in Step 1, and we have applied [27, Lemma 1] allowing to change the reference 
projector in the trace. In ([70]) . we have used the fact that Qa — Qa satisfies 
-Pa < Qa~ Qa < I ~ Pa- Wc havc also uscd that trpo {(P" - \~^)Qa} 
is uniformly bounded since it converges to trpo {{D^ ~ D * \ ■ |^^)(P — P^)}, as 
proved in the first step. By and dTO)) . we infer that 

tr{|I?" -9*1- \-'\{Qa - Qa)'} + ^D{pq^,pqJ < C 

for some uniform constant C . Using (|59p and Lemma [T71 we eventually deduce 
that {Qa) is bounded in &2{^a)- Now, we can write 

S^^AQc) > ti-.poJiV'' ~i>*\-\-')Qa} + 0{a) (71) 
> inf trpo {(X>° - I • |-i)Q} + 0(a) (72) 

Q£Qa 
tipo {Q)=N 
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since the exchange term is 0{a) and the direct term is > 0. Next we have 
hm inf trpo{(2?0-PH-r')Q} = /'(^) 

a^O Q£Qa 

trpO (Q)=N 

defined in (|52p . It suffices to compute the above infimum by a formula similar 
to (155]) and to use the convergence of the spectrum of — P * | • to the one 
of D'' - * I • This shows that lim inf ^^Idf (A^) > I'^iN) which ends 
the proof of 



Step 4: Existence of a minimizer for a small enough. We now assume that 
D satisfies the assumptions (a) and (6) of Theorem [51 Since I?° — tD * \ ■ 
has no eigenvalue which crosses when t E [0; 1], one classically deduces that 
qo = trpo (P - pO) = 0. Hence, r{N) = J2f=i ^hen TV > 0. On the other 
hand, I°{N) = \N\ for aU N E Z. This shows that {H{) is satisfied for the 
noninteracting linear model obtained in the limit a — s- 0. To prove that {H[) 
holds for a small enough is not difficult. We just have to prove that only 
finitely many strict inequalities have to be checked in {H[), and then to apply 
(f24|) . Unfortunately, we cannot use Lemma [3] since the lower bound of ([30|l 
diverges when v — vja. Instead, we prove the following 

Lemma 18. We assume that D £ C is such that ker(L'*' — D * \ ■ |^^) = {0}. 
Then there exists < ao < 4/7r and positive constants Ki, K2 such that, for any 
a € [0;ao], E'^^°'{N) > ki\N\—K2. Therefore, there exists a positive constant Kq 
independent of a such that E''/"{N) < E^/^ iN - K) + E"{K) for all \K\>Ka. 

Proof. We argue like in Step 3. Let Q e Qa be a state such that S^^^piQ) < 
E^/°'[N) + 77 and tr^o (Q) = N for some fixed ry > 0. Then by 

(1 - a7r/4) tr(|I?0|Q2) < g^(Q)\N\ + ^D[D, P) + r; 

la 



and therefore 



a /■ /■ \Q{x,y)\ 



2 



TT 



SyjRe \x-y\ 4(l-a7r/4) 
We obtain 

fBDF(Q)>tl>o(I?"Q)-I?(z7,PQ) 

TT 



dxdy<— —{goiO)a\N\+D{iy,P)/2 + r)a). 



-igo{0)a\N\+Dii?,i?)/2 + rja). 



4(1 - an/A) 
For a small enough, we have 

= tTj,o{{V° - ly ^ \ ■ \-^)Q) 

- trp^ {(V" \-'){Q + Q„)) - tr^o ((p" r')Qc.) 

> {k/2)\ trpo (g + Q,)| - trpo ((p" ly ^ \ ■ \-^)Q^) 

> k\N\/2 - k\ trpo (g„)|/2 - trpo i{V° - * | ' r')Qa) 
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(recall that Qa is defined in ([Sn]) ). We deduce that 



for some uniform constant K2. Recall that 5o(0) implicitly depends on a, and 
that it converges to 1 as a ^ 0, see 001 [29]. □ □ 

Once we know that there exists a minimizer Qa of E'^/°'{N), it is an easy 
adaptation of the previous arguments to prove that Qa takes the form ()25p and 
behaves as stated. This ends the proof of Theorem [2l □ 



6 Proof of Theorem [3] 

Step 1: Scaling properties and the spectrum of when c ^ 1. To avoid 
any confusion, we shall use the following notation in the proof: we denote by 
^a.c a(^) the infimum of the BDF energy in Qa(-^), depending on the coupling 
constant a, the speed of light c and the ultraviolet cut-off A. We are then 
interested in the limit of E'( ^ ^ ^ as c — > oo. Since most of our previous results 
are expressed in terms of the coupling constant a — 1/c, we shall often use in 
this proof the following obvious scaling property 

El,^^^,{N)=c'E-(j^^^,^^{N), (73) 
with Vcix) = c^^v{x / c). More precisely, we introduce the following operator 

^ (C/cV)(a;) =c-3/2^(a;/c). ^ > 

Then for any state Q e QcAoiN), Qc UcQU* belongs to Qao(^), and 
^i.c.cAgiQ) ~ ^^^i/c 1 Aai^'^Q^*)- '^'^ avoid both any confusion and any com- 
plicated notation, we shall always denote by and the free mean-field 
operator and the free projector when a = 1, A = cAq and the speed of light is 
c. For the other equivalent units where a = 1/c, we use the following notation: 

V°:=Ucrlu: and 2?°^ , ^ 

It will be implicit below that T)'^, T'^ and indeed all depend on c and 
Aq. We similarly write 

V^ip) = gi{\p\)a ■ LUp + 5o(b|)/3, f>^{p) = gi(b|)a • ujp + ~go{\p\)P 

and notice that go and gi are uniformly bounded with respect to a = 1/c and 
satisfy go{x) = c^go{x/c) and gi{x) = c^gi{x/c). For c large enough, we are 
able to identify the essential spectrum threshold of as stated in the following 

Lemma 19. Assume that c is large enough, then there exist k,k' > depending 
only on Aq such that 



VgoiOy + Kc-^W < < ffo(O) + (1 + '^'/c)^ (75) 

for any p E -6(0, cAq). In particular, mina{\'D'^\) — .go(0). 
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Proof. It is known that p i-^ T^^{p) is smooth hence p i— > 50 (bl) and 

P '-^ 5i(bl)/bl are also smooth. By [1T1[23], there exist two continuous functions 
ho and hi, uniformly bounded on 5(0, Aq) with respect to c, such that 

5o(bl)-.9o(0) + ^/io(bl), ffi(b|) = H(5'i(0) + ^/ii(b|)), (76) 
c c 

where 5o(0) = 1 + 0{l/c) and 5^(0) = 1 + 0(l/c). We infer from ^ that 

5o(bl) = c^3o(0) + ^/io(bl/c), gi(bl) = c\pU{0) + ^/ii(|p|/c)). (77) 
c 

Therefore there exist two positive constants kq and ki such that (7o(|p|)^ ^ 
c^9o{0)^ " KocbP and 5i(|p|)^ > Kic^^p, hence 

5o(bl)' > c^goiOf + c2|pp(Ki - /«o/c) > go(0)2 + Kc^^p 

for c large enough. Similarly, we notice that, by [29, Theorem 2.2] 

V9o{\p\r+9i{\p\r-9o{0) < .9o(bl)yTTHW-5o(0) 
which ends the proof of jTS]) since go(0)/c2 ^ -^^q^ ^ ^ 0{l/c). □ □ 



Step 2: Upper bound. Let us start by proving the upper bound 

limsup{£;r,,,,A„(^) - NgoiO)} < E^p{N). (78) 

C — >CXJ 

Let c„ — > oo be a sequence which realizes the limsup in (j78p . Let = 
(?Ai, V'Tv) be a minimizer of the Hartree-Fock energy [UHS], belonging to 
i7^(]R^, C^)^. We introduce the following subspace of Sj^ 

W„ Span = 1,.. .,7V}, ip, = 

By ([77|) . {VlfiTV^ifij) = Sij +0(l/c„) and we can thus choose an orthonormal 
basis ((/s", ip^) of W„ which satisfies \\ip" — Villj/i ^ as c„ ^ oo. We then 
take 7" = \ f?){'P?\ as a trial state. Using ([75)) . we infer 

AT 

trpo(2?V)-iVffo(0) = E ((^ffodPl)' + .91(IPI)' - .90(0)) p^, p^) 

1=1 

1=1 

This allows to prove the upper bound (|78p . 

Step 3: Lower bound: construction of an approximate solution. The main part 
of the proof will now consist in showing the lower bound 

liminf{K,.,A„(^) - Ng^m > E^p{N), (79) 
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which wiU end the proof of (j29p . To this end, we consider a sequence c„ — > oo 
which reahzes the hm inf in (j79p . For any c„ , we shall need a state Q„ which is 
not only an approximate minimizer of E'[ c„ c„Ao(-^)' ^^^'^ approximate 
solution of the self-consistent equation. Such a state will be obtained by a 
general perturbation result due to Borwein and Preis [8] (see also |23]), and 
which we state in the simplified Hilbert case as follows: 

Theorem 4 (A smooth variational perturbation principle [SII23])- Let Ai be a 

closed subset of a Hilbert space Ti., and F : A4 i-^ (—00; +00] be a lower semi- 
continuous function that is bounded from below and not identical to +00. For 
all e > and all u G M such that F{u) < inf^vi F + e^, there exist v A4 and 
w £ Conv(A^) such that 

1. F{v) < infM F + e^ 

2. \\u - v\\-^ < yfl, \v - < ^/l 

3. + e ||w - w||^ = min + e ||z - , zeAlj. 

We apply this result by taking F — £'( c„Ao' ~ ©2(-^c„Ao)i £ — and 
X := {p-po e 62(i3c„Ao) I P = P2 ^P*, trpo(P-pO) =iv} 

C Qc„A„(A^). 

We recall that the BDF functional £1 c„Ao continuous on Qc„ Ao(-^) for 

the &i (i5c„Ao) topology, hence on M. for the ©2(-^c„Ao) topology. One has 
-^i,c„ c„Ao(-^) = ii^^A^ ^i',c„,c„Ao Lieb's variational principle. Proposition [T] 
Notice also that ConvA^ = Qc„Ao {^)- Applying TheoremlH we therefore obtain 
an orthogonal projector P" on i3c„Ao and a state P" G Qc„Ao(^) such that 
Qn pn _-pO g _yvi c Qc„Ao (-^) minimizes the following perturbed functional 
QcM^ £1^^^^^;^^{Q) + "^iv{{Q - R^f} on M and satisfies 

^r,c„,c„Ao(Q") < Sr.c„.c„Ao(A^) + HQ" " ^" I e. (i5.„A„ ) < <^n' ■ 

Noticing that 

tr{(Q" ~ P")2} = 2trpo {Q"(1/2 -Vl- P")} + tr((P")2), 

since tr((g")2) ^ tr((Q")++ - (Q")"") = 2trpo {Q"(l/2 - -p°)}, it is then 
an easy adaptation of Proposition [2] to prove that Q" satisfies the following 
equation, for some /i" G M, 

g» + 7,0 ^ pn ^ + |.(l/2 _ pO _ pn)^ ^ (80) 

We then introduce the approximate vacuum solution 

Qvac := X(-oo;0) (^^Q" + |-(l/2 - - P")) - P° 

and the approximate electronic solution 7" := Q" — (5"ac- 
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step 3: Estimates on the approximate vacuum solution (3"ac- To apply previous 
results, we now introduce Q" :— Uc^Q^U*^, Qvac •= f^c„Qvacf^c„ where Uc^ is 
the scaling operator defined above in ((7^ . One has 

Qvac = X(-oo;0) + C-^Pq^ - ^.J H ' " ^n'^^ 

2 ~ 2 



r.4 



i?" + 4(1/2 - ^" )) -P^!. (81) 



Notice the obvious property V_ = Xi-oo-a) 

(v" + ^(1/2 - )) for n large 
enough, and that D{vc^,Vc„) = Cn^D{v,v). Since Q" satisfies 



i^l/c„,l,Ao(^) < fl/c,.,l,Ao(Q") < Sl/c„,l,Ao(^) + C~ 
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(Q")n>i is bounded uniformly in &^ {^Aq), for n large enough by Lemma[TJ 
In particular, pg,, is uniformly bounded in C for n large enough. We deduce 

from the equation ((8T|) satisfied by Qvac the results of [27l [28] that 
tr ((1 + \Vm:,,fY^ + ||pq„^^ = 0(l/c„). 

Using now 

tr((Q:;ac)') = tr((Q::ac)') and D{pQ._, pQ._) = c„ D(pq„^^ , Pq„^J, 
we eventually obtain 

IQvaclle.(i5. A„)=0(l/c„) and tr (|V|(Q:^,J2)1/2 _^ ^ ^(^-1/2)^ 

(82) 



Step 4: Non-relativistic limit of the approximate electronic solution 7" and 
proof of the lower bound ((75)) . By ((5^ . we have ||P" — ^ ^ < 1 for n 

large enough, and therefore that the vacuum has a vanishing charge [571 Lemma 
2]: tr-po (Qvac) = 0- Since by construction the full state Q" has a total charge 
N > 0, this means that necessarily /x" in l|80p is a positive real constant, and 
that the perturbed mean- field operator 'Dqr, + 2(1/2 — V''! — i?")/ has at least 
N positive eigenvalues. The operator 7" is then the projector on the N first 
positive eigenstates: we can write 7" — J2iLi where each ip"^ is a 

solution of the following equation 




(A^Dili being the N first positive eigenvalues of Vq^ + ^(1/2 -V^- -R"). In 
order to prove the lower bound (17^ . we shall now show that $" — ((^", ip^) 
converges to a solution of the Hartree-Fock equations. To this end, we use ideas 
from Esteban and Sere [20]: we prove that ($") is bounded in H^{R^,C^)'^ and 
that each pf stays away from the essential spectrum of Pgn + ^ (1/2 — 7^° — i?") 

as n grows. We then apply a result of Lions [33] . 
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Lemma 20. There exists a constant e > depending only on N , v and Kq, 
such that 

Vi = l,...,iV, limsup(^7 -5o(0)) < -e < 0. (84) 

n — 'oo 

Proof. First, we notice that for any i = 1, iV, /i" is at most the A^th eigenvalue 
of the following operator, with the self-interaction removed: 



\x - y\ 



with = E,^. - i?{^.y) - E,^.'^"(^)v',"(y)*, 



\x-v\ 



-(1/2-7'° -i?"). 



We estimate by means of the min-max characterization of the eigenvalues 
in the gap which was proposed by Dolbeault, Esteban and Sere in [1^]. By 
the continuation principle of I13j, one can prove that the assumptions of |131 
Theorem 1] are satisfied for c„ large enough, and therefore the first N eigenvalues 
of I?" are given by the formula 



^ik{V^)^ inf sup {V^ip,^), 



k = 1,...,7V. 



(85) 



dim V—k 



Ik II 



We then argue like in [TUl Lemma 4.5] to estimate ^n{T^^) > /^"- We choose 
an A^-dimensional vector subspace W of _ff^(M^,R), of smooth radial functions 
with compact support in the Fourier domain, and introduce 






V y 



It is clear that V^Vr is an iV-dimensional vector space for c„ large enough, 
uniformly in i? > 1. We then use ff_ © r%VR =S)^_ + Vr to estimate ^J-nCD^) 
by Formula ((55)l . Let he ip £ Sj^ + Vr such that if — p- + x with ip- g Sj^ and 
X e Vr. We compute 

First, we use Kato's inequality to obtain 



\x - y\ 



■dx dy 



< 2 IIP. 



< 



^ IlLi (|V|(/3-,(p- 

7r(Z + A^- 1) 



The same argument with Hardy's inequality leads to 

p?{x)p-{y)x{x) 



\x - y\ 



-dx dy 



<2{Z + N-1)\\p4\^. IVxIL. 

2 



/ II l|2 , 
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for some constant ki > independent of c„. Similarly, we write 



L2 



7,0p0|l/2^ 



L2 



Then, we notice that 

p0|p0|l/2 2 

By Lemma [TOl 



v/.9o(bl)2+.9i(bl)2-.9o(bl)<2.C2b|^ 



(86) 



for some constant K2 > depending only on Aq and for c„ large enough, which 
proves that 

1/T.O \| ^ /i-nOl \l/2||^ II ^ I , K3 II ||2 

\{T^^-,X)\ < K2{\T>'\ip-,ip-) |Vxii2 < ^ + 1x1^2 . 

We now estimate the term {T>iX,x)- First we use (f77)) to obtain 



(2?°x,x) = (so(bl)x,x) < ^7o(0) llxlli^ + ^ llx 



llL2 ■ 



Then, we write 



// 



pr(a^)x(;/)^ 
k - y\ 



■dxdy < (A^- 1) 



xivf 
\y\ 



■dy 



where we have used that x is a radial function. On the other hand, being 
fixed in L^, one has 

JjRe \x-y\ Jr3 |y| 

Eventually, we estimate the term involving /C„ . We use 



\L0 



and a similar inequality for Q"g^^{x , y) / \x — y| to prove that for some constant 
Ke > 0, |/C„p < (k6)^(1 — A)/c„, and therefore |/C„| < KeCn^^^\^l — A. Using 
the same method as above to estimate the term {ICnip, 'p) and the fact that 
7f(x, y)/\x — y\ defines a nonnegative operator, we therefore obtain the bound 



(2?r^,^) < (-l/2 + O(c-i/2))(|I?0|^_,^_) + 2||^_| 



2 

L2 



R 



iixiii2 



Finally, by means of < I^pAl^ II^-X|L2 < ^^s ||</3-|Il2 || X II l2 /(c„i?), we 

conclude that there exists a constant e > depending only on Aq, v, N and 
the chosen space W, such that for R large enough, < (50(0) — e + 

0(cn^^^)) ||<i3||i2 ■ This ends the proof of Lemma 



□ 



□ 
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Lemma 21. Each $" = ((^f , ...,>f%) is bounded in H^{R^,C-^)^ as n-^ oo. 

Proof. We adapt arguments from [20] . Using the self-consistent equation ((83|) 
and estimates similar to those of the proof of Lemma [201 one can prove that 
there exists a constant £ > (independent of c„, but depending on N, v and 
Aq) such that 

tr(|P°|V) < iV5o(0)2 +nr((-A)7") +^c2 tr((-A)7")i/2. 
Using now (|75p. we infer that 

tr(|I?"|V)>^ffo(0)^ + Aj4tr((-A)7"), 
which shows that X)^! llV^sf ||^2 = tr((-A)7") is bounded. □ □ 

The sequences ((/j") being bounded in i/^(R'^, C*) for all i = 1, iV, we can 
now rewrite the self-consistent equation ([55)) as 



I ■ I F-y| / 

where lim„^oo lk"llH-i(K3 c*) ~ ^ ^^"^ ()77p . We now apply the method of 
[20] to conclude that ((^", ip^) converges towards ((/?i, (/Sjv) with = , 
■0 G i/^(M'^,C^)^ which is a solution of the Hartree-Fock equations, and that 

lim {£r,c„,c„Ao(7") - ^5o(0)} = ^SfItv-)- 

By the estimates of the proof of Lemma and ([5^ . we deduce that 



n^co y JJ^e \x-y\ J 

i.e. that the interaction between the vacuum and the rest of the system vanishes. 
Since c„.c„Ao(*3vac) ^ by the stability of the free vacuum, we finally obtain 
the lower bound ([79]) . This ends the proof of ([29]) . 

Step 5: Conclusion. For c large enough > 50 (0) hence by Lemma [31 

(1 - 4/(c7r))go(0)k| - ^D{i., v) < E''{N) < 5o(0)k| 

This implies that, for c large enough, E^^N) > E^iN ~ K) + E"{K) for any 
K <0 or K > N. On the other hand, it is well-known [JTJ [33] that 

E^^iN) < mm{E^piN - K) + E^p{K), K ^ 1, N}. 

This proves that {H[) holds for c large enough, by (|29p . Thus there exists a 
minimizer Qc for E'^ {N) . It then suffices to apply again the analysis of Steps 
3-4 to show that Qc satisfies ((MJ), and obtain the stated convergence of the 
electronic orbitals towards a minimizer of the Hartree-Fock energy. This ends 
the proof of Theorem [31 □ 
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A Proof of Lemma [T] 



We write as usual Q = Q^^ + Q^^ + Q^^ + Q ^, where by assumption 
(3++,(3 e &i{Sja) and Q+~,Q^+ G &2{^a)- First, using ^ we see that 
||PQ|li2 < '^A HQ II 62(^3 a) some constant ka = 0{A^^^). Therefore, we now 
only estimate D{pQ, pq) in terms of ||Q||]^.po . The diagonal terms are treated 
as follows 

rrr^ dk<\\pQ++\\^^ \k\ dk 

B(o,A) rl Jb{o,A} 

< (2.)3(47rA) ||PQ++||1,(«3) < (2^)^(4^A) ||Q++fe,(^,) ■ 

For the off-diagonal terms, we use ideas from 27, p. 540-547]. Let C be a 
function in C fl i^A- We compute 



|(PQ+-,C)| = |tr(Q+-C)l = |tr(Q+-7'°Cn)l 

/-po I 

II62(^ia) 



^l|Q^1le2(.A) ll^-Cni 



Let us now fix some ao < A/tt. Similarly to [571 Lemma 12], it can be proved 
that there exists a positive constant k independent of a (but depending on A 
and ao) such that trc4 (7^° (?3)'P!jL(g)) < K\p- q\'^. Therefore 

IKCKLa^,) < iVCii. - 4^«' IICIIc' and \\pq,- I < n" ||Q+- L^c^a) 
for some constant k" depending only on A and ao- □ 



B On the Structure of the Variational Set 

In this section, we consider an infinite-dimensional Hilbert space and a refer- 
ence orthogonal projector lionS) such that H and 1 — n are both of infinite rank. 
We introduce ^+ := (1 — 11)^3 and := WS). First we prove a useful reduction 
for projectors belonging to H -I- 6p(i3) (i.e. to n -I- 62(^3) by [27l Lemma 2]). 
This decomposition is valid in a more general setting (for any Fredholm pair of 
projections (P, H) [2^) but for the sake of simplicity, we restrict ourselves to the 
Hilbert-Schmidt case needed in the article. Then, we deduce a general structure 
result for the variational set 

Q:={Qe6?(io) |Q* = 0, -n<Q<i-n}. (89) 

Theorem 5 (Projections and BDF-states in the Fock space). Let P he an 
orthogonal projector on such that P — H € &2{S^). Denote by (/i, .... /at) G 
{Sj+)^ an orthonormal basis of Ei = ker(P — 11 — 1) = ker(n) n ker(l — P) 
and by {gi, ...,gM) G (^-)*^ an orthonormal basis of E-i = ker(P — 11 + 1) = 
ker(l — n) n ker(P). Then there exist an orthonormal basis (wi)i>i C io+ of 
(Ei)^ in Sj+, an orthonormal basis (ui)i>i C S}- of (-B_i)^ in Sj^, and a 
sequence {Xi)i>i G ^2(1^^) such that 

P - E \fn){fn\ + E + + (90) 

n=l i—1 
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1 P I \/ I , \vi ~ Kui){vi - XiUi\ 

1 - = 2^ l5m)(5m| + 2^ • (91) 

m— 1 i—1 * 

The Bogoliubov-Dirac-Fock state 1271 associated with P in the Fock space T built 
on the electron-positron decomposition = i5+ © is then given by 

N M 



np = kY[a*ifn)Y[b*(.9m)eMAa*b*)^o (92) 

n— 1 m— 1 

N M 

= fc n «s(/«) n ^s(3™) n (1 + A.as(«.)6S(«o) (93) 



n— 1 m— 1 'i>l 

w/iere A := X]i>i and k = n.t>i(l + -^f)"^''^- 

Formula is classical and can be found in different forms in [5^[5^[FT1[5^ 
(see also [BJ Theorem 2.2]). Since 

M 

n= X! \9n,){gm\+^\Ui){Uil 
m—l i>l 

we obtain 

N M -2 

n— 1 m—l i> 1 ^ 

+ Etta^(I"^)("»i + I"')("'I)- (95) 

i>l * 

The terms in (|94l) form the diagonal part of P — H, which is trace-class [27l 
Lemma 2]. The last term (|95l) is the off-diagonal term which is only Hilbert- 
Schmidt a priori. Note we obtain from this formula that trn(P — 11) = — M 
is an integer [271 Lemma 2] . The formula of P — H can also be written as 

N M 

P = n + El/n)(/n|- E \9m){gm\+Q{A) 
n—1 m—l 

where A '}2i>i M'Vi){ui\ and 

A* A AA* ^ I 1 

Q{A) = — 7 rir- + A — - H — -A*. 

^ 1 + A*A l + AA* 1 + A*A 1 + A*A 

Therefore 

{p I p = p* = p2, p - n e 62(io), IIP - n|| < 1} 

= {Il + Q{A), Aee2iSj+,Sj-)}. (96) 

Proof of Theorem [3 We only sketch the proof which is an easy adaptation of 
ideas in [21 |36l |59l ESI jSH |5T] . Let U he a. unitary transformation such that 

p = c/nc/-i. We introduce U++ = (1 - n)[/(i - n), [/+_ = (1 - n)[/n, 

[/_+ = nC/(l - n) and = HUH. it can be verified that U+- and U-+ 
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are Hilbert-Schmidt operators, and that Ei — kerf/^^ and E^i = kerC/*_. 
The operator : ker(t/__)^ Ran(?7__) = kci-{Ul_)^ ~ E^^ possesses an 
inverse Uz]_ well-defined and bounded on E^^. Following [59i Equation (10.84)], 

we introduce the Hilbert-Schmidt operator A :— UZ- : E^^ E^. It can 

be proved that Ran(P) = Ei ffi^ (1 + A){E^^) which means that {Ei)^ is the 
graph of A. Writing A = X)i>i where {\i)i>i G £2(1^'^), one obtains 

that {Ei)^ = (1 + ^){E^i) — spanjwi -I- \iVi, i > 1}, and therefore 



The same argument applies to (j?T|) . The proof that, in the Fock space based on 
the decomposition — ^3+ (B the dressed vacuum Op is given by formula 
([92]) is let to the reader. Recall that Op is characterized by the normalization 
constraint ||0p||^ = 1 and the relations ap(/)Op = bp{f)ilp = for all f G Sj, 



whereap(/) = ao((l-P)/)-H6S((l-P)/)and6p(/)=a*(P/)+6o(P/). □ □ 



We can now clarify the structure of the variational set Q defined in ((89|) . 

Theorem 6 (Structure of the Variational Set). The set Q coincides with the 
set containing all the operators of the form 



1. D e 62(^5) is such that kei D D kerH and kei D* D ker(l - U); 

2. Ud = exp(P> 

5. 7 S ©1(^3) is a self-adjoint and trace-class operator such that [7,11] = 
and, denoting 7 ~ 11711 and 7++ = (1 — 11)7(1 — 11), then —11 < 7 < 
and < 7++ < 1 - H. 

Proof. Notice first that any Q of the form (|97p belongs to Q. Indeed UdiU^d G 
61 (^) and Ud^U-d — 11 is a difference of two orthogonal projectors which is 
in &2{^) since D g 62(^0) and therefore belongs to &^{^) by [23 Lemma 2]. 
The constraint —11 < Q < 1 — 11 is obviously satisfied. We now prove that any 
Q £ Q can be written as in ([57)) . 

Lemma 22. For any Q Cz Q, there exists an orthogonal projector P and a 
trace-class operator 7' such that [P, 7'] — and 



Moreover, P and 7' can be chosen such that trn(P — H) = 0. 

Proof of Lemma[2R Let be Q G Q. Since Q is compact, the essential spectrum 
of Q + n is {0, 1}. We write Q-f-H = X(i/2;i](Q + n)-|-7' and show that tr \y\ < 
00, which will prove ([55)1 . We can find an orthonormal basis {(pn)n>i U (V'n)n>i 
of such that 




Q = UD{n + j)U-D-n 



(97) 



where 



Q = p-n + 7'. 



(98) 



Q + n = ^ r„|(/9„)((^„| -f ^(1 - Sn)\tpn){'ipn\ 



n>l n>l 
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where r„ G [0; 1/2], s„ G [0; 1/2) and lim„^oo fn = lim„^oo s„ = 0. Computing 

tr(g++ - Q—) = J2 rn 1(1 - n)^„f + - ^n) liny^ni' 

n>l n>l 

+ ^(1 - s„) 1(1 - n)^p„f + ^ s„ inv-nll' 

n>l n>l 

which is finite for Q E G^i^), we get 

||n(y5„f < oo, ^ 1(1 - n)V'„|P < oo and ^(r„ + s„) < oo. 

n>l n>l n>l 

Then 7' = J2n>i^n\<Pn) {fn\ " 1 s« I ^n) (V'n | bsfongs to the trace class 
Si (55). The above decomposition can be easily modified in order to ensure 
that trn(P — H) = 0. It suffices to add eigenstates of Q + 11 in [0, 1/2] to 
X(i/2;i](Q + n) if trn(-P — n) < 0, or to remove eigenstates of Q + 11 in (1/2, 1] 
tox(i/2;i](Q + n) if trn(P-n) > 0. □ □ 

Let us now write g = P - n + 7' with 7' e 61(53), [7', P] = and trn(P - 
n) — 0. We apply Theorem O to P and find an operator A = X]i>i & 
&2{Sj) and two orthonormal systems (/i,.-.,/Ar) G (-^+)^ and (gi,...,gjv) G 
i^-)^ such that ([50]) and ([M]) hold. Let us then introduce e_fe = gk and 
Gfe = /fe for k — 1, A^, e-N-j = and CAr+j = vj for j > 1. We also define 
the Bogoliubov angles [11] as 9k = 7r/2 for fc = 1, N and O^^j ~ arccos(l + 
A2)-i/2 for j > 1. Then P = UdHU-d with D = Efc>i ^fc|e-fc)(efc|. Thus, 
introducing 7 = U-DYUn, we obtain the result Q = C/_D(n+7){7_£)— 11. □ □ 
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